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ABSTRACT 

This paper presents the Hubble Ultra Deep Field (HUDF), a one milhon second exposure of 
an 11 square minute-of-arc region in the southern sky with the Advanced Camera for Surveys on 
the Hubble Space Telescope using Director's Discretionary Time. The exposure time was divided 
among four filters, F435W (B435), F606W (Veoe), F775W (1775), and F850LP (zgso), to give 
approximately uniform limiting magnitudes mAB ^ 29 for point sources. The image contains at 
least 10,000 objects presented here as a catalog, the vast majority of which are galaxies. Visual 
inspection of the images shows few if any galaxies at redshifts greater than ~ 4 that resemble 
present day spiral or elliptical galaxies. The image reinforces the conclusion from the original 
Hubble Deep Field that galaxies evolved strongly during the first few billion years in the infancy 
of the universe. Using the Lyman break dropout method to derive samples of galaxies at redshifts 
between 4 and 7, it is possible to study the apparent evolution of the galaxy luminosity function 
and number density. Examination of the catalog for dropout sources yields 504 B435-dropouts, 
204 Veoe-dropouts, and 54 i775-dropouts. The i775-dropouts are most likely galaxies at redshifts 
between 6 and 7. Using these samples that are at different redshifts but derived from the same 
data, we find no evidence for a change in the characteristic luminosity of galaxies but some evi- 
dence for a decrease in their number densities between redshifts of 4 and 7. Assessing the factors 
needed to derive the luminosity function from the data suggests there is considerable uncertainty 
in parameters from samples discovered with different instruments and derived using independent 
assumptions about the source populations. This assessment calls into question some of the strong 
conclusions of recently published work on distant galaxies. The ultraviolet luminosity density of 
these samples is dominated by galaxies fainter than the characteristic luminosity, and the HUDF 
reveals considerably more luminosity than shallower surveys. The apparent ultraviolet luminos- 
ity density of galaxies appears to decrease from redshifts of a few to redshifts greater than 6, 
although this decrease may be the result of faint-end incompleteness in the most distant samples. 
The highest redshift samples show that star formation was already vigorous at the earliest epochs 
that galaxies have been observed, less than one billion years after the Big Bang. 

Subject headings: astronomical data bases: miscellaneous — cosmology: early universe — galaxies: 
evolution — galaxies: high-redshift 
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1. INTRODUCTION 

A primary motivation for deep exposures of the 
sky has been to detect the most distant objects al- 
lowed by the observing technology. Over the last 
ten years, the use of ground-based telescopes com- 
bined with the Hubble Space Telescope produced 
large samples of galaxies at redshifts as high as 5 
to study early structure formation and the assem- 
bly of stars into present-day galaxies (Steidel et al. 
1996a,b, 1999; Ellis 1998; Giavalisco 2002). These 
programs successfully revealed the distant popula- 
tions recognized for several decades as important 
for understanding how the present-day universe 
came to be (Eggen et al. 1962; Partridge and Pee- 
bles 1967a, b; Tinsley 1972a, b). Because of the 
complications arising from star formation, gas dy- 
namics, and feedback into the early intergalactic 
medium, theoretical predictions about the earliest 
galaxies are challenging, and the subject has been 
driven mainly by observations. 

Even though it has been possible to detect 
galaxies at redshifts above one, it has been diffi- 
cult to determine the redshifts and thus distances 
to objects from images only, where large samples 
may be rapidly assembled. Early workers rec- 
ognized that Lja radiation should be especially 
prominent around the first generation of galax- 
ies, despite some uncertainty about the amount 
of scattering and absorption, and there should be 
a strong edge or break in the rest frame UV spec- 
tra at 912 A owing to absorption by hydrogen 
internal to the galaxies and in the intergalactic 
medium (e.g. Partridge 1974; Davis & Wilkinson 
1974; Koo & Kron 1980.) Subsequently, Steidel 
and Hamilton (Steidel & Hamilton 1992; Steidel 
1996a, b) developed search techniques to exploit 
the Lyman edge using broad band colors to find 
galaxies with a paucity of short-wavelength flux, 
the so-called "dropout" galaxies. This technique 
has proven most productive in discovering large 
samples of high redshift galaxies in multi-band im- 
ages. There are now samples of several thousand 
galaxies at redshifts between about 2 and 5 (Stei- 
del et al. 1999, hereafter SAGDP99; Steidel et al. 
2003; Giavalisco et al. 2004.) 

When it became evident that the most dis- 
tant galaxies were characterized by compact high- 
surface brightness features (Driver et al. 1995), 
the Hubble Space Telescope took a prominent role 



in the study of young galaxies. An important ad- 
vance came from the Hubble Deep Field (HDF; 
Williams et al. 1996), a four-band, 0.5 million 
second exposure with the Wide Field Planetary 
Camera 2. This seminal program using 150 orbits 
of Director's Discretionary time on Hubble uncov- 
ered a large number of sources at redshifts above 1 
that would have been difficult to discover from the 
ground. The HDF revealed a population of small, 
irregular galaxies that often appeared in pairs or 
small groups. Much of the light from these objects 
was high surface brightness — owing to high rates 
of star formation — but concentrated, requiring the 
resolution of Hubble to identify them as distant 
galaxies as opposed to red stars, say. Extension 
of the deep field approach to the southern hemi- 
sphere (Williams et al. 2000) confirmed the main 
conclusions of the HDF but also showed the limi- 
tations of a pencil beam survey in drawing broad 
conclusions about distant populations; cosmic er- 
ror within small fields can be substantial. 

Several advances since the HDF suggested that 
even deeper observations could reveal important 
aspects of the way that galaxies were created. 
The standard cosmology holds that the atoms in 
the universe were neutral following recombination 
at a redshift, z ^ 1100, until somewhere around 
z ~ 6 — 10, at which time they were reionized by 
stars and black holes. The first observation of this 
epoch came with detection of the Gunn-Peterson 
hydrogen edge in the spectra of distant quasars 
discovered in the Sloan Digital Sky Survey (Becker 
et al. 2001, Fan et al. 2002), putting the redshift 
of complete reionization around 6. The WMAP 
experiment made an indirect determination of a 
reionization era that started as early as redshift 
11 (Kogut et al. 2003; Spergel et al. 2003; Spergel 
et al. 2006). Cold Dark Matter (CDM) models 
with a cosmological constant had some constraints 
that were not in accord with such an early epoch 
of reionization (Frenk et al. 1985), although most 
of these models have sufficient freedom to accom- 
modate even the most discrepant data. Precise 
determination of the reionization history of the 
universe remains one of the important goals of ob- 
servational astronomy. 

The luminosity function inferred from the HDF 
suggested that searching a wider area to less depth 
than the HDF would be efficient at picking up 
large populations of high redshift galaxies. An 
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important advance since the HDF was the Great 
Observatories Origins Deep Survey: GOODS (Gi- 
avalisco et al. 2004.) GOODS used the Advanced 
Camera for Surveys (ACS; Ford et al. 2003) 
on Hubble to image an area thirty times larger 
but 1 magnitude shallower than the HDF. The 
GOODS sample contains more than 60,000 galax- 
ies with photometric magnitudes in four bands, 
B435 (F445W), V6O6 (F606W), irrs (F775W), and 
Z850 (F850LP), and sufficient resolution to study 
structures as small as 1 kpc at redshifts approach- 
ing 6. This sample is excellent for statistical stud- 
ies of bright galaxies at high redshifts. 

Deep fields have an advantage over shallow 
fields to study the faint end of the luminosity func- 
tion and for increasing sample sizes when the slope 
of the luminosity function is large near the limit- 
ing magnitude of the survey. For a steep slope, the 
sample size will increase faster by investing addi- 
tional observing time in more exposure on a single 
field rather than covering more area. The luminos- 
ity function at high redshifts is imprecise, but the 
current evidence indicates that it is consistent with 
a Schechter function (Schechter 1976) with a char- 
acteristic luminosity, L*, somewhat brighter than 
the local value and a faint end slope steep enough 
to warrant investment in a deep field (SAGDP99, 
Gabasch 2004a,b.) 

The redshift at which the limiting magnitude of 
GOODS makes a deep field preferable to a wide 
field can be estimated using a standard Schechter 
function. A deep field becomes preferable at a 
redshift greater than 5, where the GOODS limit- 
ing magnitude is ^ L^, depending on the exact 
assumptions about how evolves with redshift. 
The upper limit to the redshifts of the objects in a 
deep optical survey is when the Lyman edge goes 
beyond the longest wavelength filter. A practical 
limit for the ACS is when intergalactic absorp- 
tion shortward of Lyo, shifts through the zgso fil- 
ter, z ^ 7.4. A deep field should produce samples 
of objects in the range 5 < z < 7 that are fainter 
than those found in GOODS and other wide sur- 
veys and allow a good characterization of the lu- 
minosity function in the early universe. 

It is most important to characterize the lu- 
minosity function below to see the transition 
from exponential to power-law form, to measure 
the slope, and to assess the total luminosity of 
faint galaxies. The GOODS survey was limited 



to studying galaxies at the bright end of the lu- 
minosity function for redshifts greater than about 
5. Even for lower redshifts, a deep field is use- 
ful to observe galaxies fainter than the character- 
istic brightness, providing important information 
about samples in the intermediate redshift ranges 
2 < z < 5 where much of the early star forma- 
tion in the universe took place. There is a strong 
degeneracy between derivations of object density 
and characteristic luminosity unless the luminos- 
ity function is well characterized below . Such a 
degeneracy hampers the interpretation of shallow 
surveys even with large samples. 

As shown in the next section, it is possible to 
reach well below out to redshifts near 7 with 
ACS on Hubble. This capability motivated a deep 
ACS field. 

The appearance of high redshift galaxies in the 
HDF and shallower surveys indicates substantial 
evolution in size and structure between early times 
and today. This evolution was already known at 
the time of the HDF, and subsequent observations 
tend to confirm the conclusion that the galaxy 
populations look markedly different at high red- 
shift compared to the present time. But the ap- 
parent morphology of high redshift galaxies is af- 
fected strongly by the loss of low-surface bright- 
ness features owing to cosmological dimming. An 
important way to test whether the loss of these 
features significantly distorts our perception of 
galaxies at high redshift is to make deeper ob- 
servations of the sample at intermediate redshifts. 
Thus, an ultra deep field can provide an impor- 
tant complement to the pioneering observations of 
the HDF, GOODS, and ground-based surveys by 
searching for low-surface brightness components of 
faint galaxies. 

It was evident that a deep field with the new 
capabilities of Hubble following the installation of 
the ACS could address several important issues in 
early galaxy formation. In addition to augment- 
ing the samples of galaxies at redshifts greater 
than 2, there was also the tantalizing possibility 
of pushing back the observational boundaries to 
redshifts greater than 6 to reach the reionization 
epoch. With these motivations in mind and fol- 
lowing the same philosophy pioneered by Robert 
Williams for the HDF, we held a series of meet- 
ings asking for advice on the scientific importance 
of another deep field and then assembled a Scien- 
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tific Advisory Committee with a wide range of ex- 
pertise to recommend specific parameters for the 
survey: choice of field, choice of filters, and depth 
needed for a meaningful advance. 

As with the original HDF, our purpose was to 
provide a public database using Director's Discre- 
tionary Time on the Hubble Space Telescope for 
community use. This paper emphasizes the pa- 
rameters of the database rather than the subse- 
quent analysis, but it provides a first-order anal- 
ysis of the data to assess changes in the galaxy 
populations from the highest observable redshifts 
until the present. 

Thus, we assembled a team at the Space Tele- 
scope Science Institute to create the deepest 
visual-band image of the universe to date and put 
the observations in the public domain for commu- 
nity analysis. Like the HDF, this is a multi-color, 
pencil beam survey in a single ACS field. We call 
the resulting multi-color image the Hubble Ultra 
Deep Field (HUDF). 

2. Observations 
2.1. Field selection 

The field choice derived from a desire to min- 
imize the celestial foreground radiation, maxi- 
mize the accessibility to other astronomical ob- 
servatories, maximize the overlap with extant or 
planned deep observations at x-ray, infrared, and 
radio wavelengths, and maximize the observing 
efficiency of Hubble. The original HDF was lo- 
cated in Hubble's continuous viewing zone (CVZ) 
to allow uninterrupted observations over a long 
period. The background light in CVZ orbits is 
often bright when observing in the part of the or- 
bit grazing the bright earth limb. The HDF over- 
came this limitation by taking images in the ul- 
traviolet filter, F300W, during the bright periods, 
because WFPC2 images in this filter are detec- 
tor noise limited and relatively unaffected by in- 
creased background. The Wide Field Camera of 
ACS is not sensitive at ultraviolet wavelengths, 
and the enhanced background of the bright CVZ 
orbits would seriously limit their usefulness. We, 
therefore, decided not to require that the target 
field be located in the CVZ, since it would not 
enhance the efficiency of the observations. 

Zodiacal dust within approximately 30° of the 
ecliptic plane is bright for Hubble; it was desirable 



to locate the field as far from the ecliptic as pos- 
sible. Declinations north of 35° are inaccessible 
from all major southern hemisphere observatories, 
particularly the planned Atacama Large Millime- 
ter Array (ALMA), designed to be an important 
tool for observations of distant galaxies. Declina- 
tions south of -40° are inaccessible from Hawaii 
and all observatories northward. At the outset, 
we concentrated on fields between -40° and +35° 
declination and more than 35° from the ecliptic 
plane. 

Within this declination range, there are a few 
places with very low Galactic dust and substan- 
tial investments of observing time from other pro- 
grams. The most prominent was the Chandra 
Deep Field South (CDF-S), a large (15^x150 field 
located in the direction 3^30"" -28°. This field 
has very low Galactic cirrus emission and atomic 
hydrogen column density (Schlegel, Finkbeiner, 
& Davis 1998), it passes through the zenith at 
the major observatories in Chile (the VLT, CTIO, 
Gemini South, Magellan, and ALMA), and it is ac- 
cessible from as far north as the VLA site in New 
Mexico. Furthermore, the CDF-S already has a 
substantial investment in deep x-ray observations 
with Chandra and XMM, and there are existing 
ACS observations through the Great Observato- 
ries Origins Deep Survey (GOODS, Giavalisco et 
al. 2004) allowing some useful comparisons for the 
HUDF. There are also deep infrared observations 
with the Spitzer Space Telescope (Dickinson et al. 
2004) 

CDF-S is larger than a single ACS field; sev- 
eral additional criteria guided the exact choice 
of pointing within it. The x-ray sensitivity with 
Chandra varies across the field, and it was de- 
sirable for the HUDF to coincide with a region 
of good x-ray sensitivity. There are several in- 
teresting objects identified through GOODS that 
deep observations would be most useful for, specif- 
ically a galaxy at redshift 5.8 and an old super- 
nova. We centered the field such that the high red- 
shift galaxy and old supernova were both covered, 
and the x-ray sensitivity was also very good. This 
choice produced a field centered on: RA (J2000) 
= 3^32"^39% Dec (J2000) = -27°4r29'M. Table 1 
lists the major characteristics of this field. 
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2.2. Filters 

The filter choice was identical to that chosen by 
the GOODS team. This choice provides enough 
color information for rudimentary classification of 
objects and enough wavelength coverage to search 
for the highest redshift galaxies. It also makes pos- 
sible an easy comparison of samples derived from 
both surveys. To detect objects at the highest 
possible redshifts, the observations needed to in- 
clude the longest wavelength filter, F850LP (zgso), 
a band that was sufficiently insensitive in WFPC2 
(F814W) to limit its use for the HDF. The adja- 
cent F775W {1775) filter gives minimal overlap but 
contiguous wavelength coverage. Together, these 
two bands provide excellent sensitivity to the high- 
est redshift objects detectable with ACS, the 1775- 
dropout sources, and are mandatory to search for 
objects at redshifts approaching 7. 

Four bands are desirable to provide crude spec- 
troscopic analysis of the objects. Since the long- 
wavelength observations are background-limited, 
additional sensitivity could be gained by adding 
images at shorter wavelengths without loss of 
signal-to-noise ratio. The Veoe filter, F606W, is 
immediately adjacent to 1775, broad, and provides 
excellent sensitivity to all objects at redshifts less 
than about 4. 

The HDF incorporated an ultraviolet filter use- 
ful for identifying dropout sources at redshifts near 
3. The ACS wide field camera is optimized for 
red wavelengths, making ultraviolet observations 
relatively insensitive. Since a primary goal was to 
identify samples at higher redshifts than the HDF, 
we chose the bluest filter to be the B435 band, 
F435W, immediately adjacent to Veoe band. Fol- 
lowing the conventions of the GOODS team, the 
bands are hereafter called B435, Veoe, irrs, and 

Z850- 

There are two other advantages to using the 
same filter set as GOODS. The overlap between 
the GOODS CDF-S field and the HUDF makes it 
possible to compare objects directly in both sur- 
veys to calibrate completeness estimates for the 
shallower survey. Furthermore, science analyses 
can be carried out on the two data sets using iden- 
tical methods and minimizing systematic differ- 
ences. 

Figure 1 plots the total detection efficiency for 
the four bands used for the observations. The fig- 



ure includes the spectrum of a model galaxy at a 
redshift of 5.8 for comparison. Figure 1 shows that 
this very high redshift galaxy produces a sharp 
drop in flux between zgso and 1775 with no flux 
at all in the shorter wavelength bands. The rela- 
tive detection efficiency also indicates the need for 
longer exposures in the longest wavelength bands 
to detect objects whose spectra are either flat or 
blue, typical of star forming galaxies at high red- 
shift. 
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Fig. 1. — The filter transmission curves for the 
four bands used in this survey (colored lines) along 
with the spectrum of a model star forming galaxy 
at z = 5.8 (black). The spectrum is from the 
model of Bruzual & Chariot (2003) for a 0.4 so- 
lar metallicity galaxy undergoing continuous star 
formation for 100 Myr using the intergalactic ab- 
sorption curves of Madau (1995). 

2.3. Depth 

There are several ways to estimate the depth 
needed to address the different problems described 
in the introduction. However, we caution at the 
outset that so little is known about objects at 
the highest redshifts of interest that a conserva- 
tive approach would require depths well beyond 
those possible even with Hubble. We recognized 
at the outset that our goal was to obtain as deep 
an observation as possible with the amount of 
discretionary time available to the Director, and 
the depth would be constrained by the available 
pool. The resulting sensitivity is, nevertheless, 
well-suited to make progress on each issue. 

We adopt the concordance cosmology, = 
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0.73, Qm = 0.27, h = 0.71, throughout this paper 
for ah analyses. 

The primary goal was to detect a statistically 
significant sample of galaxies at redshifts between 
5 and 7. We set this goal at about 100 objects. 
Two additional goals are to study the luminos- 
ity function of high redshift galaxies at the faint 
end down to ~ O.IL*, and to observe low surface 
brightness features in galaxies that are missed in 
shallower surveys such as GOODS and the HDF, 
both of which will be aided by the faintest lim- 
iting magnitude that can be achieved by Hubble. 
The following paragraphs estimate what might be 
achieved with a limiting AB magnitude of ~ 29"^, 
say. 

Assuming galaxy luminosity distributions are 
described by a Schechter function, it is possible 
to estimate the number of galaxies accessible to 
observation in any volume element of the uni- 
verse. SAGDP99 determined a characteristic ap- 
parent magnitude at redshift 3, m^{3)^ of 24.5 (B- 
band) in the rest frame ultraviolet, correspond- 
ing to = —20.8 corrected to the assumed cos- 
mology; the local value is —20.2 (Schechter 1976.) 
Assuming = -20.2, = 0.016 Mpc"^, and 
a = —1.6, the number of galaxies detected as a 
function of redshift per unit redshift in a single 
field of AQ sr (11 arcmin^ for ACS) is: 
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dV 
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.L'^e-^dL, (1) 



where ^ is the co-moving volume, and the frac- 
tion, X, of detectable with limiting magnitude. 



is: 



X(^Z) = RLi^l()iM*+kbiz)-miim,)/2.5 ^ (^2) 

(1 + z) 

with the luminosity distance, I^l, in units of 10 pc 
and ki){z) the k-correction at redshift, z. These 
equations may be evaluated numerically; the re- 
sults are given in Table 2 for an assumed magni- 
tude limit between 28 and 29. We also calculated 
the expected number of objects appearing in the 
Veoe, 1775, and zgso bands assuming all objects had 
the spectrum of a source undergoing continuous 
star formation for 100 Myr (Bruzual and Chariot 
2003) and accounting for intergalactic hydrogen 
absorption (Madau 1995) to illustrate where the 
objects drop out of each filter. 



The numbers in Table 2 are certainly much 
larger than a real survey would see. We have not 
attempted to correct for many observational ef- 
fects that would preclude galaxy detection in a 
survey or for the expected mixture of source sizes, 
types, colors, etc. Furthermore, the density, 0*, 
of distant populations appears to be several times 
smaller than the local value used in these calcu- 
lations (SAGDP99.) The point of this estimate is 
to show that the expected number of sources is so 
large that even a few percent detection probability 
would yield statistically useful samples of galaxies. 

It is evident from Table 2 that even a detection 
limit of 28^ in zgso would easily satisfy the goal 
of detecting ~ 100 objects above redshift 5. It 
would also reach below at redshifts approach- 
ing 7. There are several large uncertainties that 
could change these numbers in either direction: if 
continues to increase beyond redshift 3 (e.g. 
Gabasch et al. 2004a, b), there will be more faint 
objects to increase the counts and move further 
into the faint end of the luminosity function; how- 
ever, if (j)^ decreases owing either to a smaller num- 
ber density of galaxies in the early universe or cos- 
mic error, the number of detected galaxies will 
decrease. The two effects offset one another for 
number counts and introduce a difficulty of inter- 
pretation without a well-characterized luminosity 
function. 

A Z850 band limiting magnitude of 28 requires of 
order one hundred orbits of dedicated Hubble ob- 
servations. The maximum available was 400, and 
those needed to be divided among the four bands 
to give adequate spectral information. We chose 
to allocate 144 orbits each to the zgso and 1775 
bands with estimated limiting magnitudes of 28.7 
and 29.2, respectively. The remaining 112 orbits 
were split equally between B435 and Veoe, whose 
estimated limiting magnitudes were then 29.1 and 
29.3. With these sensitivities, the HUDF would, 
therefore, move firmly into the range needed to as- 
semble samples of high redshift galaxies and might 
even see substantial evolution in the properties of 
galaxies when compared with later epochs. 

2.4. Schedule and field orientation 

The observations were scheduled for two peri- 
ods during which the roll angles could be con- 
trolled to produce a nearly square image. Schedul- 
ing 400 orbits at the same pointing and with con- 
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strained orientations required the use of four roll 
angles in all: 40°, 44°, 310°, and 314° for the po- 
sition angle of the +U3 axis on the sky to increase 
the target visibility and facilitate scheduling. Ta- 
ble 3 lists the schedule of observations in orbits for 
each filter. 

Each observation or visit consisted of two orbits 
with two exposures per orbit. The exposure time 
was typically 1200 seconds but in a few cases the 
exposures had to be shortened to 850 seconds. The 
total exposure time is just under 1 million seconds. 

2.5. Small Scale Pointing: Dithers 

In addition to the large rotations between dif- 
ferent phases of the observations that were re- 
quired for scheduling, smaller shifts in telescope 
pointing were applied to different observations at 
the same position angle. These small changes in 
pointing between exposures, referred to as dithers, 
were introduced on two scales. 

First, small-scale dithers were applied to each 
of the four exposures within a two orbit visit. 
This dither pattern improved the sampling of 
the final image by introducing half-pixel off- 
sets. The ACS/WFC detector critically samples 
the point spread function (PSF) in the reddest 
bands but significantly undersamples the PSF at 
shorter wavelengths. Such undersampling leads to 
loss of spatial information and aliasing artifacts. 
The introduction of sub-pixel dithering improves 
the sampling and allows the reconstruction of a 
higher-resolution final image and a reduction of 
artifacts. 

In the case of ACS/WFC, half-pixel dithers or 
small integer numbers of pixels plus a half-pixel 
in both X and Y directions provide adequate sam- 
pling. The integer pixel components of the dither- 
ing were chosen to create the most compact dither 
pattern that ensured that a bad row or column 
could not overlap in the combined image because 
the pointings were always at least 1.5 pixels away 
from others in both X and Y. This final four-point 
dither, suggested by Stefano Casertano, is given 
in Table 4. The most compact pattern was chosen 
because the exact sub-pixel shifts will be different 
far from the center of the detectors owing to the 
very large non-linear component of the ACS/WFC 
optical distortion. The dither pattern minimizes 
this effect and ensures good sampling across the 



field. 

Additional dithers of approximately 3 and 6 
arcseconds in length in the direction perpendic- 
ular to the gap between the two ACS/WFC chips 
were introduced between visits. These offsets en- 
sure that the regions of sky falling in the gap be- 
tween the two ACS/WFC chips had at least two 
thirds of the exposure of the rest of the field and 
hence minimized the lack of uniformity of the final 
exposure map. 

3. Data analysis 

3.1. Basic data reduction 

Each of the ACS/WFC exposures was pro- 
cessed through the standard ACS data pipeline, 
CALACS. The first step removed the bias level, 
subtracted the dark current, corrected the fiat field 
and gain variations, eliminated known bad pixels, 
and calculated the photometric zeropoint. 

To achieve optimal calibration, several refer- 
ence files were created specifically for these ob- 
servations: improved dark current correction files 
(hyperdarks), improved flatfields, and bad pixel 
files. 

The hyperdarks were created using all the dark- 
current frames from the 6-month period encom- 
passing the HUDF observations. These files pro- 
vide higher signal-to-noise ratios than the typi- 
cal dark reference files that are subtracted during 
standard calibration and provide a more accurate 
representation of the overall dark current struc- 
ture appropriate to the HUDF exposures. 

New flatfield images for each filter were pro- 
duced by applying a flatfield technique that cor- 
rects only low spatial-frequency variations based 
on stellar photometry of 47 Tucanae. These new 
flatfield images (L-flats) produced a more uni- 
form sky level across the images than the stan- 
dard pipeline products. After re-calibrating the 
data with these L-flats, the images had residual 
flux of order 2% of the sky level that we ascribe 
to scattered light from bright sources outside the 
field of view. We produced images of these residu- 
als from the re-calibrated images and subsequently 
removed the scattered light by the following pro- 
cedure: 

1. All exposures from the pipeline were com- 
bined to create an image of the field. 
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2. All astronomical objects were identified in 
this image to create a mask that eliminated 
all pixels with detectable light from an ob- 
ject; this process also provided improved cos- 
mic ray masks. 

3. A median image was created from the origi- 
nal calibrated exposures using the combined 
object masks and cosmic ray plus bad pixel 
masks to exclude all pixels affected by celes- 
tial sources, cosmic rays or bad pixels. 

4. The median image was convolved with a 100 
pixel-wide smoothing function correspond- 
ing to the spatial frequency of the flat-field 
residuals. 

5. The calibrated flies were divided by this 
smoothed median image to produce another 
image of the sky. 

6. This sky image was scaled and subtracted 
from each individual calibrated image to re- 
move the scattered light. 

There is a bad-pixel flle speciflc to the HUDF 
data, containing a number of additional bad 
columns and other defects that were not present 
in the standard data quality arrays when the data 
were taken. The bad pixel flle for each band was 
created after correcting the images with the im- 
proved sky values in the previous step, then sub- 
tracting an image made by combining all expo- 
sures registered to the same pointing. The result- 
ing images made it easy to identify bad pixels, 
cosmic rays, and satellite tracks that were signif- 
icantly above the noise level. A difference image 
was then created for each exposure to determine 
the root-mean-square variations of each pixel in- 
cluding Poisson noise from any objects. All 800 
images were subsequently combined using both 
median and averaging to produce an image that 
showed the median or average deviation of each 
pixel. Additional bad pixels were identifled as 
those exceeding a threshold set to flve times the 
variance of the distribution ("sigma-clipping"). 
This technique was especially valuable in identi- 
fying charge traps not present in the original data 
quality flle. 

Finally, the bias level correction performed by 
CALACS did not fully remove the bias levels in the 
four amplifler quadrants but left residual offsets of 



a few tenths of a photo-electron. These residuals 
were corrected by flrst reversing the multiplica- 
tive flatfleld correction, solving iteratively for the 
residual bias differences between the quadrants, 
and removing these differences before re-applying 
the flatfleld. 

The ACS/WFC geometric distortion calibra- 
tion developed by King (2005) given for ACS by 
Anderson and King (ISR 2006-001) from large 
numbers of images of the globular cluster 47 Tu- 
canae provided a means to reflne the astromet- 
ric coordinates of the flnal images. This cali- 
bration includes fllter-dependent scale variations, 
new fourth-order correction polynomials and ad- 
ditional distortion correction images to model re- 
maining systematic effects. The latter cannot be 
modeled by polynomials of lower order but can in- 
troduce shifts of up to 0.2 ACS/WFC pixels. The 
inclusion of these corrections into the drizzle soft- 
ware (Fruchter & Hook 2002) that was used for 
the image combination has reduced the RMS ge- 
ometric correction error to signiflcantly less than 
0.1 pixels across the full fleld. 

3.2. Image combination 

The calibrated images were combined into a 
single image for each fllter by means of the Mul- 
tiDrizzle program (Koekemoer 2002). This pro- 
gram flrst performs sky subtraction on each in- 
put exposure, after which it uses the "drizzle" ap- 
proach (Fruchter & Hook 2002) with tools devel- 
oped for the HUDF to correct for the geometric 
distortion in each exposure and remove the shifts 
and rotations that were introduced by the obser- 
vational dither pattern, thereby producing a set of 
geometrically rectifled output images that are all 
registered onto a common grid. 

The MultiDrizzle program then combined all 
the rectifled images for each band to create a clean 
median image after having rejected the highest 
and lowest values at each pixel to help minimize 
the effect of cosmic rays and negative- valued pix- 
els. The clean median image was then transformed 
back to the original distorted frame of each input 
exposure to identify cosmic rays. 

The algorithm for identifying cosmic rays de- 
pends upon comparing the original input image, 
/in, with the clean image. Id, as well as the deriva- 
tive of the clean image. Aid, to reject pixels devi- 
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ating from the clean image using a technique simi- 
lar to that employed for the Hubble Deep Field but 
with parameters tuned specifically for the HUDF 
(cf. Williams et al. 1996). The algorithm identi- 
fies cosmic rays using the following criterion: 

\kn - Icl\ > SAIcl + ^^^(jln^G\Icl^Isky\ 

(3) 

where S was set to 1.5, S/N was set to 3.5, (Jm 
is the read-out noise, G is the electronic A-to-D 
conversion ratio, and I sky is the sky value mea- 
sured in the original input image. This procedure 
is essentially sigma-clipping, with the addition of 
the derivative image which helps to soften the re- 
jection and prevents pixels from being incorrectly 
identified as cosmic rays in regions of extremely 
sharp gradients, such as near bright stars or galaxy 
cores. A second iteration using a lower level of 
statistical significance for the detection threshold 
identified additional pixels to reject surrounding 
those that were rejected in the first iteration. This 
technique ensured robust rejection of cosmic rays 
while at the same time avoiding over-rejection of 
useful but noisy pixels in bright objects. 

The resulting cosmic ray masks were used as 
input to the final drizzle combination of all the 
images. Each image was weighted by the inverse 
variance of the exposure at the mean sky level as 
calculated from the noise model before combina- 
tion. This choice of weighting ensured that the in- 
verse variance of the final image is the sum of the 
inverse variances of the input images. An input 
variance image was calculated separately for each 
exposure, taking into account the fiatfield varia- 
tion, the sky level, the read noise and dark current, 
and all the bad pixel information. 

This step essentially performs a weighted sum 
of the input images and allows input pixels to be 
shrunk by a specific amount before being mapped 
onto the output plane. The final pixel scale was 
set to 0.6 of the input ACS pixels or O.OSO'Ypixel. 
Since this scale provides Nyquist-limited sampling 
of the PSF in all 4 bands, it is the optimum pixel 
size to use. It was also the same as used for the 
GOODS data products, allowing direct compar- 
isons of the data sets without transformation. 

A number of tests provided a means to opti- 
mize the choice of drizzle parameters, in particu- 
lar to examine the effect of the drizzle kernel on 
the uniformity of the output weight maps and the 



extent of the noise correlation in the final drizzled 
image. Drizzle provides the ability to shrink the 
input pixels before mapping them onto the out- 
put plane, which reduces the degree of correlated 
noise but can lead to large variations in the out- 
put weight map if there are insufficient numbers 
of images. We examined the weight map statistics 
in detail for a wide range of kernel sizes (given by 
the value of the drizzle parameter pixfrac) aim- 
ing to ensure that pixfrac values ranging from 0.3 
down to did not introduce too large a degree of 
variation in the weight map. 

Since the HUDF has such a large number of 
pointings in each band with a good sampling of 
the sub-pixel space, the weight map was uni- 
form and well behaved across the entire central 
region of the field, including the intersection of 
the chip gaps which typically had half the expo- 
sure coverage of the remainder of the field, even 
with pixfrac = 0. Therefore, the final images 
were drizzled using a point kernel (i.e. setting 
pixfrac = 0), corresponding to pure interlacing, 
so that each input pixel maps onto only a single 
output pixel. This choice has the advantage of 
minimizing the amount by which the output im- 
age has been processed (convolved with a smooth- 
ing function), thus providing the sharpest possible 
image and the least amount of correlated noise. 

3.3. Data quality assessments 

A series of tests confirmed the quality of the 
point-spread-function (PSF) in the final combined 
images. These aimed to quantify the extent to 
which the original Hubble PSF could be recovered, 
as well as verifying that the astrometric alignment 
of the full set of 800 images was sufficiently ac- 
curate not to affect the resulting PSF. A sample 
of isolated stars was identified across the entire 
HUDF image, covering a range of locations and 
brightnesses. PSF fits produced radial profiles of 
these stars as well as determining quantities such 
as their full width at half maximum (FWHM), 
enclosed fiux and Moffat function fit parameters. 
The mean FWHM of the stars was 0''.084 in B435, 
0'^079 in V6O6, 0'^081 in 1775 and 0'^089 in zgso, 
with a scatter of 1-2 milli-seconds of arc about 
these values. These values agree to within 2 sigma 
with the expected diffraction-limited values for 
Hubble, after taking into account the initial convo- 
lution due to discrete sampling by the 0'^05 ACS 
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pixels, as well as PSF smearing by the charge dif- 
fusion kernel between adjacent pixels and subse- 
quent convolution by the 0'^03 output pixel size. 
Figure 2 presents the PSF at the center and at a 
position near the edge of the image. 




0.05 0.1 0.15 0.2 0.25 

Radius (arcsec) 



Fig. 2. — The point-spread function derived from 
stars in the image shown at two positions: near 
the center of the image (solid line) and near the 
edge of the image (filled circles). 

The final image noise may be compared with 
the sky-limited predictions in two ways. The 
first is to calculate the noise in small regions 
that are free from obvious sources and compare 
with the estimated values from typical sky bright- 
ness. This was done for 50-pixel areas, equiva- 
lent to a square aperture 0" .2 on a side. The 
predicted background-limited variances are cal- 
culated from the sky counts in each filter as 

-2.5 \0g{SN v/50 {Rskyte^p + ^^read)/W)+^^. 

where SN is the signal-to-noise ratio (assumed to 
be 10), 50 is the number of pixels in the circular 
aperture, Rsky is the count rate in photo-electrons 
s~^ pix~^ from the sky, texp is the exposure time, 
N is the number of readouts in texp (102 for B435 
& Veoe, 288 for 1775 & zgso), cTread = 5 is the read 
out noise, and ZP is the zero point magnitude 
(AB.) This calculation gives the magnitude of a 
uniform disk of 50-pixel area whose fiux would be 
10 times the noise in the image across the same 
area. 



Table 5 compares the results (10<J/0'^2) with 
expected sensitivities calculated from the aver- 
age sky values, readout noises, and the zero 
points from Table 6.1 in the ACS Handbook 
(http: / / www.stsci.edu/hst/acs/documents/hand 
books/cyclel4/c06_expcalc3.html#328554). The 
magnitudes corresponding to 10 times the rms 
noise of a 50-pixel area achieved in the HUDF are 
very close to those predicted for purely zodiacal- 
light limited performance. The small differences 
between the predictions and the results are likely 
due to the use of an average sky value for the pre- 
dictions as opposed to the actual sky value in the 
direction of the HUDF at the time of the observa- 
tions. These estimates demonstrate that the ACS 
continued to gain sensitivity with the square root 
of the integration time even in exposures of more 
than 340 ksec. The observations achieved the nat- 
ural limits allowed by zodiacal emission and the 
overall transmission and quantum efficiency of the 
instrument. 

The second method mimics the technique used 
to construct the source catalogs by adding ar- 
tificial sources of known brightness and size to 
the data sets and using the SExtractor program 
(Bertin & Arnouts 1996) to recover them. These 
Monte-Carlo tests on the zgso image simulated 
the recovery of actual galaxies by using input 
sources similar to elliptical ( "bulge-like" ) and spi- 
ral ("disk-like") galaxies. The bulge-like sources in 
this simulation were oblate spheroids with a uni- 
form distribution of minor /major axes (^) from 
0.3 to 0.9 at arbitrary projection angles on the sky. 
The disk-like sources included a variety of orien- 
tations for thin disks, and the results discussed 
below represent an average over the orientations. 
Figure 3 shows the results of the recovery simula- 
tions, plotted as contours with the recovery frac- 
tions as functions of input source radius (major 
axis for disks) and total source magnitude. 

In these simulations, there was no correction for 
those cases where the artificial sources overlapped 
with actual objects in the image, in which case 
SExtractor typically did not recover the artificial 
sources separately. The effect was to decrease the 
recovery fraction as the source magnitude became 
fainter. It is most pronounced for the extended 
objects that have a higher chance of overlap. The 
effect is apparent in the contour plots giving re- 
covery fractions less than 1 even for objects much 
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Fig. 3. — Recovery fractions for bulge-like (a) and disk-like (b) sources are shown as a function of input 
magnitude and input source radius in zgso. Contour levels starting from the large magnitudes are: 0.1 
through 0.9 in steps of 0.1. The recovery fractions are uncorrected for the cases where the input sources 
overlapped objects in the image, leading to a slight decrease in recovery fraction as both the source size and 
magnitude increase. 



brighter than the limiting magnitude. 

The simulations make it straightforward to de- 
rive the limiting zgso magnitude as a function of 
input source size. The contour at which 50% of 
the galaxies are recovered is well described by 
quadratic fits for input source radii, r (seconds 
of arc) , between 0^' and V^: 

mSObuige = 29.27 - 3.98r + 1.91r^ (4) 
mSOdisk = 29.60 -6.40r + 1.99r^ (5) 

The sensitivity to average surface brightness for 
actual galaxies also varies with source size, at least 
for sources smaller than in radius, the ones of 
most interest for studying very high redshift ob- 
jects. 

Table 5 also shows the limiting magnitude at 
50% recovery point from the Monte-Carlo simu- 
lations for a point source in zgso ("50% Recov- 
ery"). The limiting sensitivity predicted this way 
is most applicable to assess the completeness of 
catalogs derived from the images. The underlying 
size distribution of sources is needed to calculate a 
robust limiting magnitude for objects in the cata- 
logs. Although galaxies with a variety of sizes up 



to were successfully recovered from these exper- 
iments, the recovered sizes were generally smaller 
than the actual artificial sources for input radii 
greater than about 0'^2, corresponding to about 
7 pixels, and the recovered magnitudes were often 
larger than the actual magnitudes as a result. 

4. Results 

Figure 4 displays a color rendition of the final 
HUDF image cropped to display an area of uni- 
form exposure. The color rendition allows one to 
distinguish red and blue galaxies easily and gives 
the viewer a good visual impression of the col- 
ors of typical galaxies in the HUDF. The URL 
archive.stsci.edu/prepds/udf/udf_hlsp.html gives 
the reader access to the complete data sets as well 
as the catalogs described below. 

4.1. Source Catalogs 

The SExtractor program (Bertin & Arnouts 
1996) produced catalogs from the final drizzled 
ACS images. We used the GOODS version of SEx- 
tractor that contains a few improvements in back- 
ground determination and segmentation; thus, the 
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Fig. 4. — The orientation is such that the left-hand edge is approximately NE and the top edge is approx- 
imately NW. The color mapping used to produce this rendition is: blue = (combination of B435 & Veoe), 
green = (combination of Veoe & 1775)7 ^ed = (combination of 1775 & zgso). See also Wherry, Blanton, &Hogg 
(2004) for a different approach to the color scheme. 
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results presented here may not be identical to 
those produced with the published version of the 
program. The inverse variance images were used 
as weights for calculating the RMS uncertainties 
associated with the flux measurements for each 
source. Two catalogs of sources were produced us- 
ing the i775 and zgso images to identify the sources 
and the "dual-image" mode in SExtractors, using 
isophotes of each source determined by the selec- 
tion band to measure the photometric brightness 
of each source in the other three bands, thereby 
producing isophotally matched magnitudes that 
can be directly compared to produce colors for 
each source. 

Catalogs were compiled independently by three 
team members using the same software with 
slightly different search parameters. The results 
were essentially the same among these indepen- 
dent trials; one version was chosen for the final 
catalog. The SExtractor parameters for this cat- 
alog were optimized for the pixel scale and PSF. 
The FWHM was set to 0''.09. Source detection 
required a minimum of 9 contiguous pixels with a 
detection threshold above 0.61 with a total of 32 
deblending sub-thresholds and a contrast parame- 
ter of 0.03. Additional details about the SExtrac- 
tor parameters are available from the electronic 
HUDF distribution. 

The i775 catalog contains 10,040 sources, at: 
http: / / archive.stsci.edu/prepds / udf/udf_hlsp.html. 
A visual inspection of the sources revealed a small 
number (< 0.1%) of spurious detections that are 
not included in the final catalog. Moreover, there 
are about 100 additional sources identified visually 
that were not picked up by SExtractor, owing to 
their proximity to brighter sources and the inade- 
quacy of the deblending algorithms. These sources 
were formally added by doing another SExtractor 
run with different deblending parameters. An 
initial list of 208 sources was produced, which 
was then reduced to a total of 100 sources af- 
ter visual inspection and rejection of sources that 
were clearly part of previously identified sources. 
Note that the SExtractor magnitudes of these 100 
sources at bands other than 1775 are suspect, owing 
to the close proximity of other bright objects. 

The Z850 catalog contains 39 sources that are 
not in the 1775 catalog. This catalog may also be 
obtained at the URL referenced above. 

While SExtractor is useful to identify and mea- 



sure the photometric properties of sources in the 
images, its use introduces some uncertainties that 
are often overlooked in subsequent analyses. Many 
objects consist of complicated structures that may 
either be lumped together as a single source or bro- 
ken up into several different sources depending on 
the choice of deblending parameter. These choices 
will then affect the distribution of source sizes, for 
example, and could lead to false conclusions about 
the nature of source sizes in different populations 
owing to the subjective criteria involved in iden- 
tifying individual sources. It could also affect the 
luminosity function by weighting either towards 
many small faint sources or fewer large brighter 
complexes combined as individual sources. 

Similarly, the photometric magnitudes and col- 
ors depend on the approach to aperture photom- 
etry. In this paper, we select sources and measure 
their colors entirely with isophotal magnitudes — 
mag-ISO in the SExtractor output — but we sub- 
sequently use the larger Kron-like aperture pho- 
tometry in the selection band — mag-AUTO — to 
measure total magnitudes. Thus, for example, 
if a source were selected for a particular sample 
in the zgso-band, the selection would be done on 
mag-ISO (z85o), the colors would be set by mag- 
ISO(band l)-mag-ISO(band 2), but the final mag- 
nitudes would be scaled to make the zgso magni- 
tude equal to mag-AUTO (zgso). 

The GOODS survey produced an image of the 
same field as the HUDF through the same filter 
set but with less exposure time. One way to check 
the completeness of catalogs produced with SEx- 
tractor is to compare the sources recovered from 
the GOODS data with the HUDF and examine 
differences in the recovered magnitudes. Figure 
5 shows the differences between the GOODS and 
HUDF catalog zgso-magnitudes as a function of 
HUDF Z850 for common sources that were detected 
by GOODS with at least 5a confidence. 

There is generally good agreement between the 
two surveys. In particular, we find no system- 
atic differences between GOODS and the UDF 
and no trend of galaxies becoming systematically 
brighter in the UDF as expected if most galax- 
ies had low-surface brightness regions lost at the 
GOODS depth and recovered by the HUDF. The 
direct comparison confirms also that the GOODS 
Z850 selected catalog is complete down to ~ 27^. 

The HUDF has the great advantage of using 
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Fig. 5. — The zgso -magnitude differences between 
GOODS and the HUDF as a function of the HUDF 
Z850- magnitudes for 1388 sources appearing in 
both catalogs. The red points are sources with 
50% flux radii less than 6 pixels (0'M8); the black 
points are sources with larger 50% flux radii. The 
dashed blue lines show the expected la scatter 
from photon, background, and detector noise only 
in a 112 pixel area corresponding to r^o =6. 

hundreds of individual images per band with ex- 
tensive dithering to reduce the observational noise 
to very low levels and provide very high signal- 
to-noise ratios on thousands of sources in the im- 
ages. Thus, the HUDF images are well suited to 
explore the uncertainties introduced into source 
samples through the choice of parameters in algo- 
rithms such as SExtractor, since the image noise 
should be negligible for the brighter objects. 

4.2. Number counts 

The number of sources grows with increas- 
ing magnitude until the limiting magnitudes are 
reached. Figure 6 plots the cumulative number of 
objects per square minute of arc as a function of 
magnitude for zgso- This plot also includes the 
number counts taken from the GOODS survey in 
the same units. 

The HUDF and GOODS number counts per 
unit area are identical in slope but offset in ab- 



Fig. 6. — Cumulative number counts of objects 
arcmin"^ as a function of magnitude for the 
HUDF (red) and the GOODS (black) surveys. 

solute value for magnitudes less than 26.5, the 
stated completeness of the GOODS data. The 
surface density of objects in the HUDF is smaller 
than in GOODS by about 10% over the range 
23 < Z850 < 26. The HUDF number counts con- 
tinue to rise until about magnitude 30 at which 
point incompleteness limits further increase. 

4.3. Lyman break galaxies 

Objects at high redshift have little to no ob- 
servable flux at wavelengths shortward of the rest- 
frame Lja line owing to the strong intergalactic 
absorption by hydrogen. When the redshift of the 
object is greater than about 3.5, the observational 
signature is a lack of flux in the B435 filter but 
otherwise detectable emission in Veoe, 17755 ^ind 
Z850- At higher redshifts, the Veoe and 1775 filter 
fluxes also diminish relative to the longer wave- 
lengths. The flux is said to "dropout" of the short 
wavelength filter (e.g. Steidel et al. 1996). 

We identified objects as dropouts according to 
the following criteria similar to those used in by 
GOODS: 

B435-dropouts (if all conditions met): 

Ba35-Vgog > (1.1 + 1406 - ^850) (6) 
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5435 -^^606 > 1.1 (7) 

^606 -^850 < 1.6 (8) 

S/N{Veo6) > 5 (9) 

S/N{ir75) > 3 (10) 

Vgoe-dropout if: 

Veoe - *775 > 1-47 + 0.89(^775 - zgso) 

or 2 (11) 

V606 -i775 > 1.2 (12) 

^775 - 2850 < 1-3 (13) 

S/N{ijj5) > 5 (14) 
S/N{Bi35) < 3 



or B435 - i775 > ^606 - *775 + 1 

(15) 

i775-dropout if: 

^775 - ^850 > 1.3 (16) 

^/A^(^85o) > 5 (17) 

^/A^(^606) < 2 

or ^606 - ^850 > 2.8 (18) 

We also used the compactness index in SEx- 
tractor, sometimes called stellarity, to reject stars 
in the images. Through trial and error by examin- 
ing different bright objects, we found that reject- 
ing sources with a stellarity greater than 0.9 (Veoe 
and Z85o) or 0.8 (1775) when the magnitudes were 
more than 2 magnitudes below the limiting magni- 
tudes effectively rejected stars without excluding 
small galaxies. 

The criteria listed in equations 6-8 were used 
to generate three lists of dropout candidates for 
the HUDF using the catalog published here. We 
inspected every source in these lists and rejected 
those that were spurious. The resulting culled lists 
contain 504 B435-dropouts, 204 Veoe-dropouts, 
and 54 ijys-dropouts over the 11 arcmin^ HUDF 
field. We applied the same criteria to the pub- 
lished GOODS catalog, culled the list visually, 
and found 1606 & 1661 B435-dropouts, 486 & 370 
Veoe-dropouts, and 117 & 125 iyrs-dropouts in the 
north (157 arcmin^) & south (161 arcmin^) fields, 
respectively. The corresponding areal densities of 
sources in the two surveys show similar ratios of 
B435:V606-i775-dropouts but are 8-10 times higher 
in the deeper HUDF image than in GOODS. This 



higher source density means we are going farther 
down the luminosity function with the additional 
depth. The added sensitivity partially compen- 
sates for the smaller HUDF area for discovering 
high-redshift objects. 

Complete tables of all dropout sources are avail- 
able at http:/ /www-int. stsci.edu/ svwb/hudf.html#man. 

In the area where GOODS overlaps the HUDF, 
there are three iyys-dropouts in GOODS and 6 in 
the HUDF, two of which are common. The one 
GOODS source not in the HUDF subset has a 
color of 1.48 in GOODS and 0.83 in the HUDF. 
The four HUDF sources not in the GOODS 
catalog have zgso magnitudes greater than 26. 
These differences can be understood because the 
GOODS i775 image is shallower than the zgso im- 
age, and noise in 1775 becomes a limiting factor for 
GOODS when selecting red 1775—2850 sources. 

Tables 6-8 list the sources in each of these sam- 
ples both for the HUDF and GOODS. The elec- 
tronic version of this paper contains the complete 
source lists; the printed tables include only the 
first 36 objects from each sample. 

Figures 7 to 9 show examples of these dropout 
sources. Each figure consists of images of 64 x 64 
pixels or 1" .92 on a side, corresponding to ap- 
proximately 12kpc in comoving coordinates. The 
dropout sources are always in the exact center of 
these images, although in some cases there are sev- 
eral objects that dropout in the small field. 

Visual inspection of the dropout sources in Fig- 
ures 7 through 9 shows a large number of irregu- 
lar structures along with many compact objects. 
High redshift galaxies appear much smaller and 
less regular than we see in the local universe, con- 
firming similar findings from other deep surveys 
with the resolution of Hubble, notably the HDF. 
The majority of dropout sources are compact, of 
order 1 kpc in extent. They often show multiple 
components and irregular structures. Few of the 
dropout sources show the regular structures of lo- 
cal spiral or elliptical galaxies. 

Lotz et al. (2006) simulated the appearance 
of four nearby galaxies as they would appear in 
the HUDF at a redshift of 4 using the ultraviolet 
images from GALEX. Figure 10 shows these sim- 
ulations alongside the 24 brightest B435-dropout 
sources from this present sample, presented as 
combinations of the Veoe- and i775-images. Few if 
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Fig. 7. — Images of the first 36 B435-dropout sources. Each image is 64 x 64 pixels (l'^92 square), corre- 
sponding to 13kpc at z = 4. The RGB colors coincide with the 1775, Veoe, and B435 bands to represent the 
source colors. The small lines denote the exact center position for the objects in the table. 
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Fig. 8. — Images of the first 36 Veoe-dropout sources. Each image is l'^92 square as in Figure 7, corresponding 
to 12kpc on a side at 2: = 5. The RGB colors coincide with the zgso, irrs, and Veoe bands. 



17 



Fig. 9. — Images of the first 36 irrs-dropout sources. Each image is l'^92 square as in Figure 7, corresponding 
to llkpc on a side at z = 6. The RGB colors coincide with the zgso, irrs, and Veoe bands. The source IDs 
in this figure are from the zsso-selected catalog and will not match with the B435- and Veoe-dropout IDs. 
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Fig. 10. — Simulated images of nearby galaxies from Lotz et al. (2006) on the left compared with images 
of the 24 brightest B435-dropout sources on the right. All figures are 3^' boxes corresponding to 21kpc at 
z = 4.1. The images on the right are a sum of Veoe- and irrs-fiux as assumed by the simulations. The 
dropout sources are denoted by the small lines in the center of each frame as in Figures 7-9. 
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any of these bright B435- dropout sources resemble 
the normal galaxies but many resemble the merger 
shown in the bottom panel. Although this simple 
comparison ignores some of the inherent problems 
of observing faint nearby galaxies with little ultra- 
violet flux (Lotz et al. 2006), the HUDF, neverthe- 
less, confirms the striking result of the HDF that 
the universe looks different at redshifts of 3 and 
above, corresponding to times when the universe 
was less than about 2 Gyr old. Elmegreen et al. 
(2005) reach a similar conclusion about galaxies in 
the HUDF by morphologically classifying all large 
galaxies in the field and comparing them with lo- 
cal distributions. 

For the iyys-dropouts, there are other source 
lists for comparison that are useful to understand 
how the combination of catalog generation and 
subsequent selection produce reproducible sam- 
ples. Bunker et al. (2004) identified 55 1775- 
dropouts in the HUDF, very similar to the 54 in 
this paper. Close examination shows that while 
there are many coincidences, the source properties 
often differ markedly owing to the different catalog 
generation processes used by the separate groups. 
The Bunker et al. sample breaks up several of our 
single sources into multiple sources and find sev- 
eral sources near bright objects that do not appear 
in our catalog owing to our conservative choice of 
deblending. Formally, only 38 of the objects in the 
current catalog are coincident with Bunker et al. 
sources to O'M. This number grows to 48 when 
the criteria for proximity is relaxed to 1" ^ showing 
that many of the sources identified here as single 
are multiple objects in the Bunker et al. cata- 
log. Of the remaining seven sources, four are close 
to bright galaxies and not separated from them 
with our choice of deblending parameter, and the 
other three do not pass our selection filter owing 
to slightly different magnitudes and uncertainties 
or observed flux in short wavelength filters. 

Yan and Windhorst (2004) find 108 1775- 
dropouts in the HUDF to a limiting magnitude 
Z85o= 30. Of these 108, 48 match STScI sources 
to O'M, a few do not appear in the STScI catalog 
owing to proximity to bright galaxies, and almost 
all the others are fainter than the limiting mag- 
nitude used in this study. Stronger deblending in 
the Yan and Windhorst sample means that their 
source sizes are typically smaller than ours, and 
several are sub-components of single sources in the 



present catalog. 

Visual inspection of all sources in both samples 
show they are real and not image artifacts. Our 
inspection of the Yan and Windhorst (2004) cata- 
log of i775 -dropout sources using similar selection 
criteria indicated that all those objects are also 
real, even though they have many objects unique 
to their catalog. The salient point is that differ- 
ent catalog parameters and different photometric 
magnitudes produced by the SExtractor software 
combine to create rather different source samples 
in this very simple case, introducing an important 
but not easily quantifiable uncertainty into subse- 
quent interpretations of the sources properties. 

4.4. Sample redshifts and volumes 

The redshift ranges sampled by these criteria 
depend on the (unknown) spectra of the sources 
and to some extent on the source magnitudes. To 
assess the impact of spectrum on selection, we 
calculated the redshift ranges over which galax- 
ies with different characteristic spectra may be se- 
lected. We adopted synthetic galaxy spectra from 
Bruzual and Chariot (2003) and convolved them 
with absorption by intergalactic hydrogen accord- 
ing to the prescription of Madau (1995) to simu- 
late galaxies at different redshifts, adding also two 
very simple spectra employing a step function at 

Figure 11 shows the different spectra in the 
rest frame (i.e. without intergalactic absorption) 
and how the different selection criteria act on the 
galaxies with different spectra as functions of red- 
shift for the three samples. Tables 9 and 10 give 
the exact redshift ranges and equivalent co-moving 
volumes per square minute of arc in the HUDF 
field without and with extinction for comparison. 
We stress that these calculations assume noiseless 
observations and assume perfect photometry, so 
they are not equivalent to calculations of the effec- 
tive volume (e.g. SAGDP99). Nevertheless, they 
are useful to illustrate how differences caused by 
spectral variations in the sources affect the sample 
selection. 

The differences in selection owing to spectral 
variations are generally modest, but they can af- 
fect the mean redshift and sample volumes espe- 
cially if reddening by dust is important. These dif- 
ferences are minimized in the highest redshift sam- 
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Fig. 11. — Upper left: Synthetic spectra of galaxies to assess dropout sample criteria. Upper right, lower 
left, lower right: Each of the criteria in equations 6-8 are plotted as left side-right side, so that they satisfy 
the inequalities when they cross zero. The quantities that increase through zero, e.g. B — V — {1.1 -\- V — z) , 
define the start of the selection range; the quantities that decrease through zero, e.g. 1.6 — (V — z), define 
the end of the selection range. The color coding for models is the same as in the upper right: lOMyr CSF 
(blue), 100 Myr CSF (green), 1000 Myr CSF (magenta). Step (black). Red step (red). The vertical "ticks" 
on the f{z) =0 line illustrate the selection redshifts {zmin and z^ax )for the step function spectrum. 
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pies. One additional effect for the iyys -dropouts is 
the importance of observational sensitivity to the 
most distant sources. The upper limit to the red- 
shift range is technically around 7.4 for very bright 
objects {L ^ L*), but in practice may be limited 
to a smaller value by the paucity of bright galaxies. 
This effect is a minor problem for the B435- and 
Veoe-dropouts where the color criteria themselves 
limit the upper end of the redshift range. 

Several papers have noted that photometric er- 
rors can scatter objects in and out of the color se- 
lection windows. This scattering will be most pro- 
nounced for selection criteria with shallow slopes 
as they approach zero in Figure 11, for exam- 
ple, f{z) = {V -i) - Min[1.47 + 0.89(z - z),2] 
corresponding to the condition 1^06 — ^775 > 
[1.47 + 0.89(z775 - Z85o) or 2] for Veoe-dropout 
limits near 2: ~ 5.3, or for the upper limit selection 
of B435 -dropouts for red spectra (cf. Tables 9 & 

10.) 

These examples do not include sources with 
strong Lyo; lines and so ignore at least one known 
population of objects in the analysis. Only very 
strong line-emission will affect the broad-band col- 
ors enough to change the selection criteria, prob- 
ably less than the variations in color seen in our 
synthetic spectra. The selection criteria in princi- 
ple allow detection of more evolved stellar popu- 
lations that are intrinsically red in redshift ranges 
outside those of interest. We tested the spectrum 
of a 1 Gyr old population at solar metallicity, also 
from Bruzual and Chariot. The B435-dropout cri- 
teria effectively excluded the selection of this pop- 
ulation at any redshift. Both the Veoe-dropout 
and i775 -dropout criteria allowed detection of some 
interlopers, but the expected contamination from 
the expected space density of such objects is small. 
In fact, the i775-dropout criteria allow detection of 
these populations over two redshift ranges, the sec- 
ond being limited by the redshift at which the age 
of the universe is l,Gyr. 

At these redshifts, the selection criteria are 
mainly sensitive to the absorption edge of the in- 
tergalactic medium for any population with de- 
tectable amounts of ultraviolet radiation. When 
there is little or no ultraviolet radiation present, as 
in a more evolved population, the criteria pick up 
the galaxies at redshifts where their intrinsic spec- 
tra show strong edges, complicating the interpre- 
tation. The dominant interlopers in these samples 



are likely to be evolved stellar populations at low 
redshifts. In addition, reddening by dust is likely 
to make some high-redshift galaxies undetectable 
in our samples, decreasing the completeness. Ex- 
tinction by dust can be the largest single uncer- 
tainty in deriving source luminosities from the lim- 
ited spectral information and can have profound 
consequences for understanding the total amount 
of light emitted by the distant populations, but 
it only affects source selection at redshifts below 
about 4 — see, for example. Figure 1 of SAGDP99. 

4.5. Size distributions 

It is possible to get a crude measure of galaxy 
sizes using the parameters produced by SExtrac- 
tor. Figure 12 plots the distributions of SExtrac- 
tor radii enclosing 50% of the flux for the three 
dropout samples in the two surveys together with 
log-normal functions for comparison. The plot- 
ted uncertainties are proportional to the square- 
root of the number of sources in each bin. Also 
shown in the figure is the size distribution for the 
1388 sources that appear in both the HUDF and 
GOODS catalogs (coordinates within O'M) mea- 
sured separately in each catalog. A log-normal dis- 
tribution characterizes the GOODS data remark- 
ably well, considering the crude measure of size 
used here. The sources in the HUDF sample ap- 
pear smaller on average for the B435- and Veoe- 
dropouts. 

The differences between the size distributions 
for the HUDF and GOODS samples indicate the 
general problem of characterizing these irregu- 
lar sources with single sizes, comparing sizes for 
different source luminosities and survey sensitivi- 
ties, and also highlights the difficulty of assessing 
changes with redshift. There is modest evidence 
for a slight increase in typical source size in the 
HUDF i775-drop sample presented here, but we 
suspect the uncertainties in deriving these sam- 
ples coupled with the difficulties of using the 50% 
flux radius to characterize them preclude any con- 
clusion that the source sizes are evolving or stay- 
ing constant. As noted in the previous section, it 
is not always straightforward to determine when a 
close grouping of several sources should be treated 
as a single, larger object. This problem resulted in 
a large mismatch between the Bunker et al. (2004) 
i775-dropout sample and the i775-dropout sample 
presented in this paper, even though the color se- 
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B-dropout distribution 



V-dropout distribution 





Fig. 12. — The size distributions of the HUDF and GOODS dropout sources and for 1388 sources in the 
GOODS catalog that also appear in the HUDF (lower right.) The size parameter is the radius (in pixels) 
enclosing 50% of the flux as returned by SExtractor. Large colored points are HUDF, small black points are 
GOODS in each plot. Uncertainties are Poisson noise of numbers per bin. The black lines are log- normal 
functions with the peak and width given as /i = log(4.6) and S = log(1.5). The colored lines have fitting 
parameters: /i = log(3.6), S = log(1.4) (blue), /i = log(3.6) and S = log(1.35) (green). 
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Fig. 13. — The assumed probability of detecting dropout sources vs. magnitude for the HUDF (left) and 
GOODS (right) in the three bands, Veoe (green), 1775 (magenta), and Z85o(red) based on the convolution of a 
log-normal size distribution with the recovery fractions as a function of size and magnitude from Figure 3. 



lection was nearly identical. 

Figures 7 — 9 show that these distant galaxies 
are often close to foreground objects that would 
be blended in images with poorer angular resolu- 
tion. A substantial fraction of these galaxies will 
be blended with foreground objects in images from 
ground-based telescopes limited to ~ f resolution 
by the Earth's atmosphere. This blending and the 
uncertainties in the models needed to correct for 
its effect introduce an additional complication in 
comparing samples derived from space-based and 
ground-based observations at faint magnitudes. 

Figure 13 shows the probability of detecting 
dropout sources as a function of source magnitude 
in the HUDF and GOODS by using the log-normal 
distribution and convolving it with the recovery 
fractions shown in Figure 3. These probabilities 
will be used subsequently to assess the expected 
number of dropout sources seen in the two sur- 
veys. 

4.6. Cosmic error 

The spatial distribution of galaxies is correlated 
over many scales, leading to fluctuations in the 
number density of galaxies. For small volumes, 
this correlation can dominate the errors in num- 
ber counts of galaxies. The cosmic error is the 
square root of the field-to-field variance due to 
large scale structure in excess of the Poisson noise 
for any sample. It is an important problem for the 
HUDF owing to the small field of view. The rela- 



tive cosmic error, a^, is the root variance corrected 
for Poisson noise, divided by the average number 
density (see Somerville et al. 2004, hereafter S04, 
equation 1). 

There are several ways to estimate this error 
for the HUDF samples. An empirical lower limit 
for the bright objects can be obtained by divid- 
ing the large GOODS area into equal-sized re- 
gions of approximately the area of the HUDF and 
looking at the variance of sources among them. 
Dividing GOODS into two 3 x 5 grids of 10.9 
square arcminute rectangles, the distribution of 
B435-dropouts among these is GOODS-N: (82, 76, 
105, 93, 134, 128, 100, 115, 115, 95, 107, 123, 116, 
110, 106), and GOODS-S: (60, 153, 140, 134, 131, 
99, 133, 108, 106, 133, 93, 95, 95, 88, 93). 

The fractional root variance of this distribu- 
tion is 19% with a statistical sampling error of 
9.6% (approximately 110 sources per cell). The 
relative cosmic error after removal of shot noise 
is then (jy = 0.16. This measurement applies 
only to the brighter B435 -dropouts detected at the 
GOODS depth, and gives a crude measure of the 
fluctuations that might be expected in the HUDF 
dropout samples. However, we note that the ac- 
tual variations among fields is as large as a factor 
of 2.5, 60 vs. 153 in the first two GOODS-S sub- 
fields, for example. 

Figure 14 plots the distribution of B435-dropout 
sources on the two GOODS fields divided into a 
3x5 grid. It is easy to see the variations among 
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Fig. 14. — The location of the B435-dropout sources are plotted on the GOODS areas. Each GOODS field is 
divided into a grid of 3 x 5 HUDF-sized areas to show the typical variance of the sample location on the sky. 
The size and color of the points depend on the Veoe magnitude as follows: large red: V606< 25.5; medium 
blue: 25.5 <V606< 27; small green: 27 <V606- 



these different HUDF-sized areas. These plots 
segment the samples into bright (red), moderate 
(blue), and faint (green) sources to allow separate 
estimations of source variance at different lumi- 
nosities. More luminous objects are more rare and 
also more clustered, leading to larger field-to-field 
variation in their number counts. 

A similar calculation for the Veoe-dropouts 
yields a relative cosmic error of 16% with a statisti- 
cal sampling error of 16%. The number of sources 
per HUDF-sized area is 32 and 25 in GOODS- 
N/S, respectively, so the field-to-field variance 
due to shot noise is comparable to the variance 
due to large scale structure. Note that the vari- 
ance between GOODS-N and GOODS-S for Veoe- 
dropouts is about 14% between two 160 square 
arcminute areas. The number of irrs-dropouts 
is too small to give a meaningful measure of the 
field- to-field variance. 

This empirical approach suffers from use of 
the brighter and potentially more clustered ob- 
jects and the use of a relatively small number of 
cells. A second approach is to use the measured 
correlation function from Lee et al. (2006) for 
GOODS B435 and Veoe- dropouts and the ana- 



lytical expression in Equation (3) of S04. For 
the faintest B435-dropout sample at limiting mag- 
nitude, Z85o= 27, the correlation scale length 
is ro = 4.14 Mpc. Assuming 7 = 1.8 implies 
cr^ (B435drop) = 0.31 for spherical cells. We cor- 
rected this result for the very elongated cell ge- 
ometry of the HUDF assuming scale independent 
bias, yielding cr^(B435drop) = 0.14. A similar es- 
timate for Veoe-dropouts using ro = 6.3, 7 = 1.8 
yields ^^(Veoedrop) = 0.49 for spherical cells and 
^^^(Veoedrop) = 0.25 for the HUDF cell geometry. 

A third approach is to use the theoretical ex- 
pectations for the clustering of dark matter in a 
hierarchical cosmology plus an estimate of the bias 
of the galaxy population in question, as described 
in S04. Here we have improved on the results pre- 
sented in S04 by correcting for the pencil-beam ge- 
ometry of the HUDF rather than using spherical 
cells. For the B435-dropouts and Veoe-dropouts, 
we can estimate the bias of the fainter HUDF 
sample by extrapolating to lower luminosities us- 
ing the halo occupation model derived by Lee et 
al. (2006) based on the brighter GOODS sam- 
ples. For the B435-dropouts, the co-moving sam- 
ple volume is 4.24 x lO^Mpc^, yielding a source 
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density of 0.011 Mpc~^ for 504 sources. From Fig- 
ure 9 of Lee et al. (2006), we estimate the bias 
at this number density to be 6 ~ 2.4. Using 
c"DM = 0.14 (for spherical cehs, as plotted in Fig- 
ure 3 of S04), we then get a y (3435 drop) = 0.34. 
After correcting for the HUDF cell geometry, this 
is reduced to (7^(8435 drop) = 0.15. Similarly, for 
Veoe-dropouts, the halo occupation model in Fig- 
ure 9 of Lee et al. (2006) and the number density 
in the HUDF gives a bias estimate of 2.7, result- 
ing in <j^(V606drop) = 0.17 with corrections to the 
HUDF geometry. 

There is no extant measure of the clustering or 
constraints on the halo occupation model for 1775- 
dropouts, but we can estimate the bias assuming 
one galaxy per halo as in S04. With a number 
density of 54/(43, 867 Mpc^) = 0.0013 at z = 6.46, 
the predicted bias is 4.7, and (7^(1775 drop) = 0.23 
for the HUDF (with the geometrical correction for 
non-spherical cells). 

The three methods give consistent estimates of 
the cosmic error for the B435- and Veoe-dropouts, 
an encouraging result considering the uncertain- 
ties. In what follows, we adopt 17% as a charac- 
teristic value for the cosmic error. 

4.7. Luminosity functions 

The distribution of apparent source bright- 
nesses results from binning the samples and com- 
paring the number density of galaxies on the sky 
at different magnitudes. It is useful to allow the 
bin sizes to vary to keep the number of objects in 
each bin similar. 

Tables 11-13 list the number of objects per 
square arcminute in different magnitude ranges 
from the samples above. The relative uncertain- 
ties are the square root of the total number of 
objects in each bin combined in quadrature with 
a fractional uncertainty of 0.17 to include at least 
minimally cosmic error. The actual uncertainties 
are certainly higher owing to selection biases es- 
pecially within the small area of the HUDF; there 
is no attempt to remove these biases here. Each 
table also lists the total surface brightness of radia- 
tion in the shortest waveband for source detection 
(e.g. Veoe for B435-dropouts.) 

Figure 15 plots these number density functions 
for the three redshift intervals of the samples. The 
GOODS N and S sources are plotted separately 



to show the variances in samples over the rela- 
tively large areas of GOODS. The previous section 
demonstrated that the empirical variance within 
an HUDF area is at least 17% for GOODS B435- 
dropouts, and this variance is likely to apply to all 
the HUDF samples. 

Figure 15 also has Schechter functions (Schech- 
ter 1976) overplotted on the data. The fits to the 
number density distributions all assume a faint 
end slope parameter, a = —1.6. The fits are done 
in two ways to illustrate the sensitivity to the pa- 
rameters and the selection biases associated with 
the sensitivity of the observations. The dashed 
lines are fits of a Schechter functional form with 
two parameters, m*, and S*, to the apparent mag- 
nitudes to derive a characteristic magnitude and 
sky density. The solid lines assume Schechter func- 
tion distributions with a characteristic absolute 
magnitude, M^^ and volume density, 0*, integrat- 
ing the function over the sample redshift interval 
weighting each redshift by the co-moving volume 
and relative contribution of the Schechter function 
for an apparent magnitude bin: 

I/* I -\- Z 

<)-^"/>-(^)"'"*^''- <^°' 

where Dl{z) is the luminosity distance at red- 
shift, z, in 10 pc units, the apparent magni- 
tude at band, 6, kb{z) is the k-correction for a 
specific spectrum, is the solid angle (arcmin^), 
a and are the usual Schechter parameters, , (j)^^ 
is the Schechter (j)^ divided by 2.51og(e), and ^ 
is the CO- moving volume at redshift, z. The red- 
shift interval, zi-zi was that appropriate to the 
Bruzual and Chariot 100 Myr continuous star for- 
mation model with 0.4 solar metallicity and with- 
out extinction (cf. Figure 11 and Table 9); dif- 
ferent model spectra mostly affected 0* owing to 
different sample volumes. The k-correction tends 
to be dominated by the position of the absorption 
edge in the relevant band in all cases, not the long 
wavelength slope of the spectrum. Indeed, a sim- 
ple Schechter function assuming all galaxies are 
at the same redshift gives almost the same best-fit 
parameters as the integral over redshift with slight 
differences in the ratios of Tj^/(j)^ among the sam- 
ples. To fix the magnitude scale, we normalized all 
galaxy spectra to make correspond to 1400 A. 
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Fig. 15. — Number counts per square minute of arc per magnitude vs. magnitude (Veoe, 1775, and zgso) for 
the dropout samples. The horizontal bars show the range of magnitudes covered for each point. The vertical 
bars indicate uncertainties that are proportional to the square root of the counts"^ in that interval plus the 
variance of 0.17. The HUDF are the light-colored points: B435-dropouts (blue), Veoe-dropouts (green), and 
i775-dropouts (red); GOODS are darker colored points that are large for GOODS-N or small for GOODS-S. 
Each sample is also compared with two Schechter function fits all with a = —1.6. The dashed lines are 
Schechter functions fitted directly to the number counts with parameters: = 26.0, = 6.6 arcmin"^ 
(B435-drop), = 26.1, = 2.4arcmin~^ (Veoa-drop), and = 26.4, = 1.7arcmin~^ (iyys-drop). The 
solid lines are in-situ Schechter functions integrated over the selection redshift intervals for 100 Myr CSF 
galaxies with parameters: = -20.3,^* = O.OOlSMpc"^ (B435-drop), = -20.4,0* = O.OOllMpc"^ 
(Veoa-drop), and M^^ = —20.5, 0* = 0.0007 Mpc~^ (iyrs-drop). The Schechter functions have been multiplied 
by the detection probabilities for the different surveys to show the anticipated dropoff in number counts; 
there have been no corrections for incompleteness applied to the data. The plot in the lower right shows the 
HST B435-dropouts compared to the G-dropouts from SAGDP99 (magenta circles) and Sawicki & Thompson 
(2006) (magenta squares) and B-band dropouts from Ouchi et al. (2004) (green squares) together with the 
fixed- Schechter fit from the upper left plot. There ai^^no completeness corrections applied to any points. 



The solid lines plotted in the figures are actu- 
ally Ps5(m5)/(m5), where Ps5(m) is the probability 
of detecting an object from sample, 5, in band, 6, 
as shown in Figure 13. It is important to note is 
that the probability of detection depends on the 
object size. Our approach assumes that the distri- 
bution of object sizes is approximately log- normal, 
as shown in the previous sections, and to use the 
same size distribution (i.e. log-normal parame- 
ters) for all three samples. This assumption only 
affects the contribution for the faintest magni- 
tude bins in the GOODS and HUDF samples, and 
the subsequent analyses use data far enough away 
from the detection limits to ensure that Vsh(j^h) is 
greater than 0.8. 

Table 14 lists the parameters for these various 
fits along with the formal normalized values. 
The fits are remarkably good considering the sim- 
plicity of the assumptions and the lack of any at- 
tempt to correct for observational biases such as 
variation of spectra within the samples, extinc- 
tion by dust, and selection bias in the presence 
of noise. This latter bias, usually modeled with 
Monte Carlo calculations, could be considerable 
at the faint end of the samples. 

To test for density evolution among different 
redshift samples, each luminosity function was fit- 
ted assuming a single Schechter function according 
to the prescription given above for integrating over 
the redshift interval and allowing only (\)^ to vary. 
The results of these fits are also given in Table 14 
in the last three columns. It was not possible to 
get acceptable fits by fixing (\)^ and allowing 
to vary. 

The conclusions of these fitting exercises can be 
summarized as follows: 

• The apparent luminosity functions may 
be fitted adequately by assuming simple 
Schechter function distributions of either the 
apparent magnitudes or intrinsic source lu- 
minosities. 

• The best fit Schechter functions have charac- 
teristic magnitudes at ~ 1400 A (rest frame) 
of -20.3 to -20.7. The local value B-band 
value is -20.2 (Schechter 1976.) 

• The characteristic volume densities, re- 
quired to fit the apparent distributions vary 
between 0.002 Mpc"^ at z - 4 to 0.0005 



Mpc ^ at z ~ 6. The local B-band value is 
0.016 Mpc-^ 

• It is possible to fit the all these luminosity 
functions adequately assuming no change in 
characteristic magnitude, M*, but requiring 
a change in the characteristic density, 
dropping by a factor of three from z ~ 4 to 
z ~ 6. 

• It not possible to fit these distributions with 
0* ~ constant and variable M^. 

• The data are consistent with a faint end 
slope a = — 1.6 and do not require variations 
of a to fit the distributions. 

• There are strong degeneracies between M*, 

and OL that allow different combinations 
to give acceptable fits, allowing for variations 
in luminosity and density among samples, 
depending on the choice of values within 
these degeneracies. 

5. Discussion 

One important goal of the HUDF was to search 
for galaxies at redshifts greater than 5 to assess 
how they have evolved in the first few billion 
years after the Big Bang. The increased sensi- 
tivity of these observations makes it possible to 
detect galaxies out to redshifts of order 7 and 
to compare samples at different redshifts looking 
for changes in characteristic luminosity, co-moving 
density, and morphology of typical objects. The 
great advantages of the HUDF are its ability to 
observe low surface brightness emission in these 
faint sources and to characterize the luminosity 
functions at magnitudes fainter than . The pri- 
mary disadvantages are the small area of the sky 
over which the observations can be made and the 
limited spectral information available for the faint 
sources, making it difficult to determine the exact 
redshifts of the objects and to eliminate many of 
the subtle biases that are unavoidable when using 
coarse filtering techniques to cull enormous sam- 
ples. 

The public release of the HUDF data stimu- 
lated considerable interest by many groups who 
used the data to study the early evolution of galax- 
ies. Indeed, one of our goals was to stimulate this 
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community interest. The subsequent research pa- 
pers cover a wide-range of techniques to analyze 
the HUDF images and discover new properties of 
galaxy populations. We encourage the reader to 
study the resulting publications to appreciate the 
richness of the data and its interpretations, includ- 
ing Bunker et al. 2004 (the first published pa- 
per); Yan and Windhorst 2004; Stiavelli, Fall, and 
Panagia 2004a, b; Stanway et al. 2004; Bouwens et 
al. 2004a, b; Thompson et al. 2004, Stanway et al. 
2005, and Straughn et al. 2006. Bouwens et al. 
(2005, hereafter BIBF05) summarizes most of the 
results of these early works and presents the most 
comprehensive analysis to date of irrs-dropouts in 
all deep ACS data from the HUDF and GOODS 
surveys as well as an extensive discussion of the 
various potential selection biases, ways to correct 
these, and results for large samples of sources from 
these data. 

Our purpose is to complement these approaches 
by adopting the simplest possible assumptions 
about the data and source samples to see if we 
can reach general conclusions that can be bounded 
by likely variations in the sample biases. In keep- 
ing with the original goals of the HUDF project, 
we confine ourselves only to the most salient re- 
sults implied by the observations. We restrict the 
analyses to samples taken from the HUDF and 
GOODS with identical selection criteria and pre- 
sumably similar observational biases. 

5.1. Evolution of galaxy luminosities 

It is of great interest to characterize the dis- 
tribution of galaxy numbers, sizes, and luminosi- 
ties at different epochs in the universe to under- 
stand how the present day distribution evolved 
from the past. One challenge to this character- 
ization is source identification itself, because the 
structures seen in observations of the early uni- 
verse are complex and do not necessarily corre- 
spond to the structures we see today. The choices 
about how sources are identified and their subse- 
quent measurements can have profound effects on 
global quantities such as the characteristic lumi- 
nosity, typical size, and number density of individ- 
ual objects. It is therefore important to separate 
as much as possible observables derived from the 
data and model-dependent assumptions about the 
early universe used to analyze these observables to 
assess the relative impacts of the assumptions on 



any interpretation of how the universe evolves. 

Source samples derived from space-based and 
ground-based observations differ in their selection 
biases and can be compared only with attention 
to these differences. It has become common prac- 
tice to correct for observational biases by deriv- 
ing an "effective volume," ]4ff, for each sample 
by simulating different selection effects in a com- 
puter and doing Monte Carlo calculations of their 
average effect on model populations with differ- 
ent luminosities, shape, and sizes at different red- 
shifts then using V^ff to correct the observed dis- 
tributions of sources to characterize the sample 
through luminosity functions (SAGDP99; Ouchi 
et al. 2004, hereafter O04; BIBF05; Sawicki & 
Thompson 2006) This method necessarily includes 
assumptions about the underlying populations of 
galaxies that are then difficult to separate from the 
observations in the final results. It is especially 
problematic when comparing samples derived us- 
ing two different observational techniques, such as 
ground-based telescopes limited by seeing and sky 
brightness with space-based telescopes, because 
the biases in each method are quite different and 
may not be equally well corrected by the effective 
volumes. 

In this paper, we present number counts per 
sohd angle, S(m), as a function of apparent mag- 
nitude as a purely observational quantity derived 
from the data (Section 4.7.) The model-dependent 
assumptions needed to fit these functions are em- 
bedded in the fitting routines explicitly. The bi- 
ases introduced by different spectral distributions 
of the underlying populations can be calculated for 
the redshift ranges allowed by the color-selection 
criteria. Although incomplete, when applied con- 
sistently to data sets gathered using the same 
method, it should be a good way to compare sam- 
ples at different redshifts without a priori knowl- 
edge of the underlying source population to dis- 
cover significant evolution of the source popula- 
tions. 

This approach suggests that the number counts 
at redshifts between 4 and 7 may be characterized 
by Schechter functions as seen in Figure 15. The 
HUDF is sufficiently deep to detect galaxies some- 
what fainter than the characteristic magnitude, 
M*, at all redshifts < 7 , a key goal of the obser- 
vations. It is therefore an excellent complement 
to shallower surveys such as GOODS that sam- 
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pie the bright end of the luminosity distribution. 
Furthermore, it is possible to fit the distributions 
well by adopting a faint end slope, a ^ —1.6, and 
a value for the characteristic absolute magnitude, 
^ —20.3, similar to the local value, although 
the wavelengths are quite different (Johnson B- 
magnitude for the local value, Mab ^ 1400 A for 
the analysis here.) 

The two salient problems in comparing the 
space-based to the ground-based data sets are that 
they sample different ends of the luminosity func- 
tion, and the comparisons are normally made after 
using different model assumptions to compute the 
effective volumes to correct observed surface den- 
sities to in situ volume densities. There are po- 
tential biases created by different model assump- 
tions in V^ff. SAGDP99 pointed out that small 
differences in the color selection criteria used for 
dropout samples can create large differences in the 
resulting absolute luminosity functions, although 
the differences are minimized when making rela- 
tive comparisons between samples observed and 
analyzed the same way. The Veoe -dropouts are 
especially sensitive to the color of the underlying 
galaxies owing to the shallow slopes of the main 
selection function, Vqoq — ^775 > 1.47 + 0.89(^775 ~ 
Z85o), near the selection boundary; see Figure 11c. 
Thus, when comparing different samples, it is im- 
portant to use the same assumptions for the anal- 
ysis of Feff. 

Figure 15 includes a comparison between the 
space-based and ground-based samples at 2: ~ 4 
from four samples: SAGDP99, O04, Sawicki and 
Thompson (2006), and the samples in this paper. 
The observed surface density data from two of the 
ground-based samples match the space-based sam- 
ples reasonably well in the regions of overlap, al- 
though the faint end of the Sawicki and Thomp- 
son data are significantly under the space-based 
and O04 counts. It is evident from this figure that 
the shape of the apparent luminosity functions be- 
tween the ground-based and space-based samples 
look different without normalization to the sample 
volumes. 

We derived rudimentary sample volumes for 
three of the samples plotted in Figure 15 adopt- 
ing the entire co-moving volumes within the color- 
selection redshifts. The absolute magnitudes are 
calculated from the apparent magnitudes using 
distance moduli and k-corrections from a Bruzual 



and Chariot continuous star formation model over 
100 Myr with 0.4 solar metallicity and no extinc- 
tion convolved with the different filter response 
functions. For the G-dropouts of SAGDP99, we 
used the spectroscopic redshift distribution for the 
volume and the ratios of Vef[ given in their Table 
4 to correct for incompleteness in the two faintest 
magnitude bins of their study. For the O04 data, 
only color-selection sample volumes were used. 
The same technique was applied to the B435- 
dropouts presented in this paper using the same 
assumptions and limiting the magnitude range to 
make completeness corrections less than 20% (cf. 
Table 9 and Figure 13.) We could not correct the 
Sawicki & Thompson (2006) data without knowl- 
edge of their completeness corrections, although 
it is evident from inspection of their figures that 
the faint end luminosity functions will not match 
the space-based counts unless the corrections for 
incompleteness are factors of 3 or more. 

Figure 16 shows the resulting B435- , B-, and G- 
dropout luminosity functions along with Schechter 
function fits to the data. The data are easy to 
fit with a single Schechter function with a char- 
acteristic luminosity and faint end slope similar 
to those derived for various samples using both 
ground and space-based observations. Figure 8 of 
SAGDP99 gives best fit parameters to the z ~ 4 
sample of = 24.97, = 0.013 Mpc"^ for an 
Einstein-de Sitter universe with h = 1. For a con- 
cordance model and k-corrections assumed in this 
paper, we calculate the equivalent = —20.9 
and (j)^ = 0.0011 Mpc~^. The values derived in 
Figure 16, = -20.7 and 0* = 0.0013 Mpc"^, 
are entirely consistent with the SAGDP99 values, 
considering the degeneracy between and 0* in 
the fits and different methods of calculating vol- 
umes between the two papers. The absolute lumi- 
nosities depend on the k-corrections applied to the 
different filter sets used for the SAGDP99, O04, 
and Hubble observations. 

This method allows us to tie the data sets to- 
gether under a single set of assumptions. However, 
the absolute values of M*, (j)^, and a depend both 
on the model assumptions and the completeness 
corrections, which we have estimated using rela- 
tively simple assumptions. It is possible to fit all 
the distributions in this paper well by assuming 
constant and varying (j)^. It is not possible to 
fit all the distributions in this paper well by as- 
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Fig. 16. — The in situ luminosity functions of the dropout sources from the HUDF, GOODS, SAGDP99, and 
O04 surveys corrected using the model-dependent Vs and distance moduli as discussed in Section 5. Color 
and shape codes are the same as in Fig. 15. Upper left: B435-dropouts. The black line is the best-fit Schechter 
function for this combined sample: = —20.7,^* = 0.0013 Mpc~^ (x^ = 0.6. The magenta line is our 
translation of the original fit from Figure 8 of SAGDP99 using = -20.9 and = 0.0011 Mpc"^ for the 
concordance cosmology and k-corrections for a 100 Myr CSF Bruzual spectrum using our model calculations. 
We were unable to include the Sawicki & Thompson (2006) data in this plot without knowledge of their 
completeness corrections. Upper right: Veoe-dropouts from the HST samples together with Veoe-dropouts 
from O04 (black squares). Fit parameters are: - 20.55,0* = 0.0009 Mpc"^ (x^ = 1.1.) Lower left: HST 
i775-dropouts. Fit parameters are: M* = —21.1,0* = 0.0005 Mpc~^ (x^ = 0.6.) 
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suming a constant density of galaxies and varying 
and a. In all cases, the average density of 
galaxies at redshifts above 6 appears smaller than 
at redshifts < 4 by a factor of order three for fixed 

Our treatment of galaxy evolution using sam- 
ples all derived with similar techniques and as- 
sumptions calls into question one of the strongest 
conclusions of BIBF05, who state "Scenarios, such 
as density evolution ((/>*), which do not include this 
evolution in M* or a are ruled out at 99.9999% 
confidence - demonstrating quite convincingly 
that galaxies at z ~ 6 are lower in luminosity 
than galaxies at z ~ 3." A simple way to test this 
statement using the samples presented here is to 
examine the ratio of luminosity density above a 
certain magnitude limit to that below the limit 
and see if this ratio changes with redshift. Ta- 
ble 16 presents the results of this test. Each 
dropout sample was divided into a bright and 
faint part with the boundary set as the apparent 
magnitude of = —20.8 galaxy at the mean 
redshift of the sample using k-corrections appro- 
priate to a Bruzual & Chariot 100 Myr continuous 
star formation model. The faint end limit was 
set 2 magnitudes above the boundary to ensure 
the inclusion of similar bright/faint magnitude 
ranges among the three samples. Table 16 gives 
the resulting surface brightness and in-situ energy 
densities as well as faint/bright ratios. It demon- 
strates that at least for these samples there is no 
evidence for an evolution of the faint/bright ratio 
within the uncertainties. 

There are certainly biases in our technique that 
could affect this conclusion. If there is less absorp- 
tion by the intergalactic medium than assumed 
by the Madau (1995) model for low redshifts (see, 
e.g., Fontanot et al. 2006), we could miss bright 
galaxies that have enough residual B435 fiux to be 
excluded as Veoe -dropouts, say, even though the 
colors might otherwise qualify: there is an extra 
filter in equation (15) to deal with this problem 
but it was determined by trial and error and may 
not be robust for large samples. BIBF05 discov- 
ered many more faint galaxies in their 1775 -dropout 
sample that will boost the faint/bright ratio in Ta- 
ble 16, an effect that would be offset by including 
the additional faint galaxies in the B435- and Veoe- 
samples by extending the faint magnitude limit 
below M -18.8. 



The biases should apply in a similar fashion 
to each of the samples presented here, however. 
The ground-based data have a different set of bi- 
ases, and it appears that even when great care 
is taken to model instrumental and observational 
effects it is difficult to calculate large corrections 
accurately. The faint end luminosity functions of 
Sawicki & Thompson (2006) in Figure 15 clearly 
show incompleteness in the last few points, even 
though their paper uses these points to establish 
a flatter faint end slope than seen in any of the 
space-based samples or in the Subaru data. Ex- 
amination of Figure 7 shows how much crowding 
exists on scales of ~ 1" ^ typical of ground-based 
seeing, and how it would be difficult to correct for 
loss of faint galaxies owing to blending with other 
objects. 

It is encouraging that Figure 16 shows a con- 
sistency between the B435-, B- and G-dropouts 
when analyzed with similar assumptions. It re- 
mains to be seen if this consistency will hold up 
under a more detailed analysis that assesses the 
impact of observational noise on the selection of 
these separate samples. Observations at longer 
wavelengths for all these samples might also be 
used to constrain the spectral variations — to assess 
the impact of dust, for example — and better reflne 
the model assumptions using the data themselves. 
Spectra of the high redshift samples derived here 
will be essential to remove some of the largest un- 
certainties. 

The actual luminosity density decreases mod- 
estly from low to high redshifts, as shown in Figure 
17 which is an update of Figure 1 from Giavalisco 
etal. (2004.) Giavalisco et al. (2004) were the flrst 
to note the difficulty of assessing effective volumes 
especially for the i775-dropout sample, where the 
sample redshift range is determined more by limit- 
ing sensitivity than pure color selection — note the 
comparison between the open and fllled circles in 
Figure 17 that use the same data but different as- 
sumptions about the effective volumes of the sam- 
ples. The HUDF greatly ameliorates this problem 
with additional sensitivity that allows detection 
of galaxies to redshifts very close to that for 
which the Lyc^ edge moves out of the zgso fllter. 
Adding the faint i775-dropout sources discovered 
by BIBF05 below our cutoff will certainly boost 
the high redshift point in this flgure, suggesting 
little if any evolution in the luminosity density be- 
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tween redshifts 6 and 3. 
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Fig. 17. — A comparison of the total energy den- 
sities from different surveys after Figure 1 of Gi- 
avalisco et al. (2004): crosses: SAGDP99, filled 
and open circles: Giavalisco et al. for two different 
T4ff , rectangles: total light from HUDF sources us- 
ing k-corrections and mean dropout distances for 
HUDF only. 

The important point is that some derived prop- 
erties of the galaxy population such as the luminos- 
ity distributions often depend more strongly on the 
assumptions used to correct biases than they do on 
the data themselves. In many cases, the correction 
factors are large, as noted by SAGDP99, mean- 
ing the absolute densities and absolute magnitudes 
can depend almost entirely on the corrections. 
Without separate consistency checks to establish 
ground truth — high resolution spectroscopy of ev- 
ery source in a field, say — it is dangerous to take 
the model results too seriously. This problem is 
exacerbated in differential comparisons between 
derived quantities from data sets taken with dif- 
ferent instruments and using assumptions about 
the underlying populations for their analysis. 

5.2. Cosmological dimming of galaxies 

For the individual objects seen in the images, a 
comparison of sources seen in both the HUDF and 
GOODS suggests that the Hubble observations 



adequately measure the total luminosity of galax- 
ies brighter than the completeness limits. Figure 5 
shows that there are few if any cases of low-surface 
brightness emission at the periphery of galaxies 
picked up in the deeper HUDF survey that are not 
seen in the shallower GOODS images. The pop- 
ulations at high redshift appear to be much more 
compact than local galaxies, making them easy to 
see with the high resolution of the Hubble images. 
The colors of the bright and faint galaxies in each 
dropout sample are similar as is the emission in 
the i775 and zgso for the B435-dropouts (cf. Fig- 
ure 18), suggesting that obscuration by dust is not 
large for the galaxies seen here. Figure 18 also in- 
dicates that the long wavelength spectra of galax- 
ies at redshifts near 4 derived using the dropout 
technique are flat, and there is not a strong indi- 
cation of changes in reddening by dust at different 
redshifts in the two lowest redshift populations. 

The apparent surface brightness of radiation 
from the entire population of high-redshift galax- 
ies is increased by the addition of the HUDF com- 
pared to the GOODS samples, typically by a fac- 
tor of about two. This higher surface brightness 
results from the contribution to the faint end of 
the luminosity functions that dominate the to- 
tal luminosity in each sample. There is no obvi- 
ous turnover in these functions (Figure 15) before 
reaching the limiting magnitudes of the HUDF 
(see also the discussion by Yan and Windhorst 
2004.) Certainly, the observed surface brightness 
is a lower limit that would be improved only by 
much deeper observations than even the HUDF. 

6. Conclusions 

The Hubble Ultra Deep Field reached limiting 
magnitudes within a few tenths of the zodiacal- 
light limited values and approximately 1.5 magni- 
tudes deeper than the next deepest observations. 
It uncovered more than 10,000 new galaxies, the 
vast majority of which are below the sensitivity 
of other surveys. Many of these sources have 
redshifts greater than three corresponding to the 
epoch less than two billion years after the Big 
Bang. 

Three samples of galaxies in different redshift 
intervals were identified in the HUDF from ob- 
jects with missing flux in the shortest wavelength 
bands, those that "dropped out" of the B435, Veoe, 
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Fig. 18. — The average spectrum for the B- and Veoe-dropout samples vs. rest-frame wavelength. All sources 
are assumed to be at the mean redshift in each sample: 4.1 (B-drops) and 5.0 (V-drops). The flux densities 
through each filter have been averaged after normalizing to the selection filter (Veoe or 1775) flux density. 
The error bars are the standard deviations of the samples. In each plot, the averages take place in three bins 
of selection magnitude. For B435-dropouts, these are: V < 26.5, 26.5 < V < 28, V > 28. The ranges are the 
same for Veoe-dropouts except the magnitude is 1775. The three sub-samples are plotted from left to right 
in the figures with decreasing symbol sizes and offset -25, 0, and 25 A from the mean rest frame wavelength. 



or i775 filters, corresponding approximately to red- 
shift ranges 3.5 < z < 4.7, 4.6 < z < 5.5, and 
5.7 < z < 7, respectively. The same technique 
was applied to the GOODS survey to get samples 
of similar objects at the bright end of the lumi- 
nosity function but over a larger region of the sky. 
Comparison of these samples showed that: 

• As found in previous observations with the 
Hubble Space Telescope, galaxies at these 
high redshifts are smaller and less symmet- 
ric in shape than galaxies at lower red- 
shifts. These results confirm the conclusion 
of the first Hubble Deep Field observations 
that typical galaxies in the early universe 
look markedly different than galaxies today, 
showing that galaxies evolved rapidly in the 
first few billion years after the Big Bang. 

• The faint sources detected by the HUDF are 
smaller on average than the brighter objects 
seen in shallower surveys such as GOODS or 
the HDF, with typical sizes (radius enclosing 
50% of the flux) about 0'M5. There is not 
strong evidence for evolution of source size 
with redshift among these samples. 

• The HUDF samples the luminosity function 



approximately 2 magnitudes below the char- 
acteristic magnitude at redshifts as large as 
6, thus fulfilling a key goal of the observa- 
tions. 

• It is possible to characterize the apparent 
luminosity functions of all samples with 
Schechter functions assuming only a vari- 
ation in source density (0*) and not lumi- 
nosity (M^) or faint-end slope {a). It is 
not possible to fit the luminosity functions 
assuming evolution of the characteristic lu- 
minosity alone. The uncertainties in de- 
riving and fitting the luminosity functions 
are large enough to accommodate luminos- 
ity evolution for comparison with theory, for 
example, but it is not required by the data. 

• The co-moving volume density of galaxies 
needed to fit these high-redshift samples is 
at least 10 times smaller than the local value 
of 0.016 Mpc~^ using the simplest assump- 
tions about source selection and the entire 
co-moving volumes allowed by color selec- 
tion. Use of an effective volume approach 
may raise these densities, and it is quite 
likely that the Lyman-break dropout tech- 
nique picks up only a fraction of the galaxies 
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at high redshift, meaning the absolute value 
of (f)^ derived here is a lower limit. 

• The shape of the bright-end luminosity func- 
tion at z ~ 4 estimated from B435-dropouts 
closely matches that of ground-based sam- 
ples using a similar dropout technique, when 
the samples are treated with similar (albeit 
uncertain) assumptions. 

• Although the luminosity density of galax- 
ies at wavelengths ~ 1400 A decreases mod- 
estly from redshifts of a few out to redshifts 
greater than 6, vigorous star formation was 
already underway when the universe was less 
then 1 billion years old. 

Our purpose in carrying out the HUDF observa- 
tions was to provide a deep dataset as a legacy to 
science for study of the early universe. The many 
publications arising from the original data release 
testify to the utility of the HUDF for this pur- 
pose. The images provide a rich hunting ground 
for those who wish to understand the nature of the 
universe at the time when stars and galaxies first 
came into existence after the Big Bang. 

These observations would not have been possi- 
ble without the help of many people at the Space 
Telescope Science Institute and in the community. 
We are especially indebted to Stefano Casertano, 
Ian Jordan, Max Mutchler, and Bill Workman. 
We profited from discussions with Mark Dickin- 
son, Richard Ellis, Michael Fall, Andy Fruchter, 
Mauro Giavalisco, Hans- Walter Rix, Dan Stark, 
Charles Steidel, Rodger Thompson, and especially 
Robert Williams. For advice on the choice of 
observational parameters as well as occasional 
discussions about the scientific directions, we 
thank the group of distinguished scientists who 
served as the scientific advisory committee to this 
project: Ron Ekers, Catherine Cesarsky, Guenther 
Hasinger, Garth Illingworth, Jeremy Mould, Matt 
Mountain, Anneila Sargent, Tom Soifer, Harvey 
Tananbaum, Robert Williams, and Rogier Wind- 
horst. We also acknowledge the many good sug- 
gestions that came out of the ACS High Latitude 
Survey Workshop in March 2001 and the Hub- 
ble Treasury workshop in Baltimore in November, 
2001. SB is grateful to the Max-Planck-Institut 
fiir Astronomic for providing research facilities and 



hosting several visits to enable much of the anal- 
ysis presented here and to a number of important 
suggestions from the MPIA staff for that analysis. 
We are also grateful to the California Institute 
for Technology for hosting one of us (SB) and 
generously sharing information about the dropout 
techniques used there. The referee was most help- 
ful in pointing out weaknesses in the arguments 
and stimulating analysis of the data. Of course, 
none of these observations could have taken place 
without the dedication of the hundreds of engi- 
neers, managers, and scientists at NASA and its 
contractors who gave us the Hubble Space Tele- 
scope as a lasting legacy to the world. 
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Table 1 
HUDF Field Characteristics 



Field Center: 



Size: 
Area: 

Position Angle: 
E(B-V): 

HI column density: 
Radio sources: 



RA(J2000) = 3 32 39.0 
Dec(J2000) = -27 47 29.1 
I = 223.6°, b = -54.4° 
Ecliptic latitude: -45.2° 
200" X 200" 
11.0 arcmin^ 
~ 42° 
0.008 

7.9 X lO^^cm-2 

two 300 mJy (365 MHz) at 

one 1 Jy (408 MHz) at 5l' 



10' 



Table 2 

Potential number of sources in HUDF 



Redshift Total ^ Veoe ^ irrs ^ 



1 26,800 23,600 23,900 13,900 

2 17,000 17,500 15,800 7,700 

3 10,700 11,500 11,000 5,300 

4 7,000 5,900 7,700 3,700 

5 4,900 930 4,800 2,600 

6 3,600 490 1,600 
7 2,700 120 



Differential numbers are rounded for clarity; mAB limits: 29, 29, 29, 28. 
Columns 3 5 represent estimates using an unreddened model spectrum 
starburst galaxy with intergalactic hydrogen absorption; 
0* = 0.016Mpc~^, = -20.2, a = -1.6. 



Table 3 
HUDF Orbits 



Date range 


PA 


B435 


Veoe 


1775 


Z850 


2003.09.24 2003.10.28 


310 


6 


8 


18 


18 


2003.09.24 2003.10.28 


314 


22 


20 


48 


56 


2003.12.04 2004.01.15 


40 


6 


8 


18 


20 


2003.12.04 2004.01.15 


44 


22 


20 


50 


50 


Total orbits 




56 


56 


144 


144 


Total ksec 




134.9 


135.3 


347.1 


346.6 



Table 4 
Dither pattern 





y" 


X (pix) 


y (pix) 


Exp # 


0.0000 


0.00000 


0.0 


0.0 


1 


0.1488 


0.08595 


3.0 


1.5 


2 


0.2232 


0.24075 


4.5 


4.5 


3 


0.0744 


0.15480 


1.5 


3.0 


4 
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Table 5 
HUDF Sensitivities 



Filter 


Exposure 


ZP 


<Sky> 


Predict 


10cr/0".2 


50% recovery 




ksec 




s""*" pix~^ 


mAB 


mAB 


mAB 


F435W 


134.9 


25.67 


0.0284 


29.09 


29.08 




F606W 


135.3 


26.49 


0.1358 


29.29 


29.45 




F775W 


347.1 


25.65 


0.0780 


29.21 


29.35 




F850LP 


346.6 


24.86 


0.0415 


28.67 


28.70 


29.4 



The ZP & <Sky> from the ACS manual. Read noise (Tread — 5. 
50% recovery magnitude is the average disk/bulge for point sources. 



Table 6 

First 36 B435-DROPOUT sources 



STScI ID 


RA (J2000) 


Dec (J2000) 


r50(") 


B 




V 




i 




z 


<Tz 


28 


53. 


16485 


-27.82724 


0. 


.089 


30.84 


0. 


.54 


29.08 


0. 


.09 


29.00 


0, 


.07 


29.50 


0.20 


44 


53. 


16797 


-27.82653 


0. 


187 


29.67 


0. 


.38 


27.92 


0. 


.08 


27.52 


0. 


.04 


27.71 


0.08 


79 


53. 


16769 


-27.82542 


0. 


115 


31.14 


0. 


.86 


28.93 


0. 


.09 


28.67 


0. 


.07 


28.69 


0.12 


91 


53. 


17054 


-27.82467 


0. 


133 


31.33 


1. 


.13 


29.02 


0. 


.21 


28.50 


0. 


.06 


28.71 


0.13 


102 


53. 


16347 


-27.82437 


0. 


116 


30.39 


0. 


.33 


29.2 


0. 


.11 


29.31 


0. 


.09 


29.28 


0.16 


148 


53. 


16454 


-27.82316 


0. 


133 


30.21 


0. 


.53 


27.53 


0. 


.04 


26.66 


0. 


.02 


26.65 


0.03 


154 


53. 


16586 


-27.82344 


0. 


167 


30.10 


0. 


.34 


28.91 


0. 


.16 


28.85 


0. 


.08 


29.60 


0.28 


230 


53. 


15400 


-27.82118 


0. 


176 


27.44 


0. 


.07 


25.67 


0. 


.01 


25.23 


0. 


.01 


25.23 


0.01 


236 


53. 


15455 


-27.82214 


0. 


.077 


31.27 


0. 


.74 


29.30 


0. 


.13 


29.03 


0. 


.07 


29.14 


0.14 


280 


53. 


15610 


-27.82122 


0. 


.121 








27.71 


0. 


.04 


27.37 


0. 


.03 


27.50 


0.05 


291 


53. 


17409 


-27.82107 


0. 


.107 


30.86 


0. 


.72 


28.44 


0. 


.07 


28.31 


0. 


.05 


28.43 


0.10 


301 


53. 


16346 


-27.82102 


0. 


.093 








29.19 


0. 


.13 


28.09 


0. 


.05 


28.22 


0.09 


307 


53. 


16511 


-27.82100 


0. 


.113 


32.12 


1. 


.74 


29.32 


0. 


.13 


29.02 


0. 


.08 


29.11 


0.15 


348 


53. 


16972 


-27.82032 


0. 


.111 


29.78 


0. 


.25 


28.32 


0. 


.06 


28.09 


0. 


.04 


28.13 


0.07 


372 


53. 


17023 


-27.81995 


0. 


.119 


30.57 


0. 


.50 


27.91 


0. 


.04 


27.87 


0. 


.03 


28.12 


0.07 


396 


53. 


17545 


-27.81957 


0. 


.140 


28.99 


0. 


.17 


27.57 


0. 


.04 


27.31 


0. 


.03 


27.24 


0.05 


401 


53. 


16197 


-27.81909 


0. 


.245 


29.23 


0. 


.39 


26.20 


0. 


.02 


25.28 


0. 


.01 


25.14 


0.01 


407 


53. 


16069 


-27.81918 


0. 


.097 


28.32 


0. 


.10 


26.11 


0. 


.01 


25.65 


0, 


.01 


25.70 


0.01 


409 


53. 


17234 


-27.81933 


0. 


.139 


28.31 


0. 


.09 


26.99 


0. 


.02 


26.83 


0, 


.02 


26.95 


0.03 


420 


53. 


15847 


-27.81947 


0. 


.094 


30.75 


0. 


.50 


28.81 


0. 


.08 


28.66 


0. 


.06 


28.67 


0.10 


421 


53. 


15600 


-27.81939 


0. 


.180 


31.44 


1. 


.26 


29.31 


0. 


.16 


28.69 


0. 


.08 


29.26 


0.22 


481 


53. 


15569 


-27.81870 


0. 


.089 


30.28 


0. 


.29 


29.06 


0. 


.09 


28.99 


0. 


.07 


29.11 


0.13 


577 


53. 


16145 


-27.81744 


0. 


.110 


30.73 


0. 


.71 


28.05 


0. 


.05 


27.09 


0. 


.02 


27.05 


0.03 


591 


53. 


17494 


-27.81738 


0. 


.108 








28.64 


0. 


.08 


28.28 


0. 


.05 


28.45 


0.09 


613 


53. 


16039 


-27.81736 


0. 


.193 


29.62 


0. 


.30 


28.04 


0. 


.07 


27.78 


0. 


.04 


27.78 


0.07 


631 


53. 


16702 


-27.81700 


0. 


.100 


29.84 


0. 


.36 


26.90 


0. 


.02 


26.43 


0. 


.01 


26.50 


0.02 


651 


53. 


14933 


-27.81696 


0. 


.084 








29.11 


0. 


.10 


28.35 


0. 


.04 


28.49 


0.08 


700 


53. 


16573 


-27.81646 


0. 


.085 


30.70 


0. 


.60 


28.27 


0. 


.05 


27.37 


0. 


.02 


27.41 


0.04 


741 


53. 


16223 


-27.81577 


0. 


.209 


27.49 


0. 


.07 


25.89 


0. 


.01 


25.45 


0. 


.01 


25.42 


0.01 


750 


53. 


17981 


-27.81593 


0. 


.120 


29.73 


0. 


.23 


28.34 


0. 


.07 


28.26 


0. 


.05 


28.37 


0.09 


751 


53. 


16381 


-27.81591 


0. 


.135 


28.59 


0. 


.11 


27.48 


0. 


.04 


27.29 


0. 


.03 


27.48 


0.05 


793 


53. 


14770 


-27.81547 


0. 


.172 


29.71 


0. 


.36 


28.03 


0. 


.07 


27.72 


0. 


.04 


27.97 


0.10 


807 


53. 


17183 


-27.81518 


0. 


.139 


30.12 


0. 


.48 


27.57 


0. 


.05 


27.15 


0. 


.03 


27.19 


0.05 


831 


53. 


14756 


-27.81495 


0. 


.107 








28.39 


0. 


.07 


27.76 


0. 


.03 


27.97 


0.07 


848 


53. 


16470 


-27.81486 


0. 


.194 


31.65 


2. 


.06 


28.05 


0. 


.11 


28.50 


0. 


.09 


28.90 


0.23 


904 


53. 


16409 


-27.81418 


0. 


.116 


30.40 


0. 


.38 


28.81 


0. 


.12 


28.33 


0. 


.04 


28.40 


0.08 
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Table 7 

First 36 Veoe-DROPOUT sources 



STScI ID 


RA (J2000) 


Dec (J2000) 


r50(") 


B 


o-B 


V 


av 


i 




z 




67 


53. 


,16471 


-^7.8z57z 


0. 


.104 


31.08 


0.69 


30.62 


0.31 


28.80 


0. 


.15 


29.00 


0.13 


123 


53, 


,16771 


TV ooono 

-z 1 .o2o)Jo 


0. 


.112 






31.21 


0.62 


28.86 


0. 


.09 


29. 


0.16 


204 


53. 


,16884 


-27.82266 


0. 


.101 


I 


I 


31.47 


0.61 


29.24 


0. 


.12 


29.61 


0.21 


384 


53. 


,15979 


-27.82000 


0. 


,066 


34.59 


12.11 


32.27 


0.99 


29.38 


0. 


,11 


29.99 


0.24 


546 


53. 


,15804 


-27.81794 


0. 


.145 






30.28 


0.35 


27.77 


0. 


.06 


27.04 


0.04 


646 


53. 


,16606 


-27.81678 


0. 


.148 






29.39 


0.14 


27.41 


0. 


.04 


27.32 


0.04 


712 


53. 


,17836 


-27.81625 


0. 


,095 


31.81 


1.87 


29.74 


0.19 


27.26 


0. 


,03 


27.11 


0.03 


748 


53. 


,17549 


-27.81601 


0. 


.138 


31.94 


1.93 


31.06 


0.60 


28.82 


0. 


.15 


28.92 


0.16 


811 


53. 


,14928 


-27.81537 


0. 


.078 






31.29 


0.47 


29.00 


0. 


.09 


29.62 


0.20 


847 


53. 


,16759 


-27.81499 


0. 


,092 






29.77 


0.15 


28.48 


0. 


.07 


28.87 


0.13 


850 


53. 


,16398 


-27.81503 


0. 


.102 






30.64 


0.25 


29.38 


0. 


.14 


29.96 


0.27 


864 


53. 


,16253 


-27.81490 


0. 


.118 


31.26 


0.80 


31.26 


0.55 


28.95 


0. 


.09 


28.93 


0.13 


947 


53. 


,16205 


-27.81389 


0. 


.082 






31.76 


0.78 


28.91 


0. 


.12 


28.97 


0.12 


971 


53. 


,16449 


-27.81373 


0. 


.064 


31.37 


0.66 


31.17 


0.38 


29.32 


0. 


.19 


29.53 


0.17 


992 


53. 


,17062 


-27.81350 


0. 


.066 


34.33 


9.19 


32.06 


0.79 


29.14 


0. 


.10 


29.71 


0.18 


1055 


53. 


,16123 


-27.81267 


0. 


.231 


32.57 


5.95 


31.93 


2.28 


28.40 


0. 


.14 


28.61 


0.21 


1126 


53. 


,17548 


-27.81208 


0. 


,097 






29.96 


0.17 


28.57 


0. 


,07 


29.11 


0.16 


1154 


53. 


,17563 


-27.81167 


0. 


.125 


31.76 


1.58 


29.50 


0.14 


28.07 


0. 


.07 


28.24 


0.09 


1156 


53. 


,16736 


-27.81171 


0. 


.351 






29.73 


0.39 


28.16 


0. 


.11 


29.00 


0.39 


1157 


53. 


,15931 


-27.81169 


0. 


.110 






30.32 


0.25 


28.39 


0. 


.06 


28.53 


0.10 


1276 


53. 


,14984 


-27.81067 


0. 


.127 






29.80 


0.19 


27.77 


0. 


.04 


27.63 


0.05 


1318 


53. 


,16177 


-27.81035 


0. 


,059 


31.34 


0.62 


30.25 


0.16 


28.94 


0. 


,08 


29.21 


0.12 


1379 


53. 


,14944 


-27.80973 


0. 


,118 


30.81 


0.65 


30.46 


0.32 


28.32 


0. 


,08 


28.47 


0.10 


1392 


53. 


,15636 


-27.80959 


0. 


,078 






30.86 


0.45 


27.63 


0. 


,03 


28.27 


0.08 


1519 


53. 


,15031 


-27.80882 


0. 


,106 






31.62 


0.67 


29.17 


0. 


,12 


29.12 


0.13 


1608 


53. 


,14378 


-27.80828 


0. 


.091 


31.62 


1.19 


33.56 


4.87 


29.01 


0. 


,10 


28.72 


0.11 


1749 


53. 


,15279 


-27.80729 


0. 


.071 






30.82 


0.27 


29.09 


0. 


,09 


29.48 


0.16 


1796 


53. 


,15324 


-27.80697 


0. 


.145 


30.95 


0.70 


30.12 


0.23 


28.21 


0. 


.06 


27.99 


0.06 


1846 


53. 


,17214 


-27.80674 


0. 


,109 






30.23 


0.23 


28.71 


0. 


.09 


29.00 


0.15 


1967 


53. 


,17512 


-27.80612 


0. 


.084 


31.78 


1.04 


30.52 


0.23 


29.29 


0. 


.15 


30.05 


0.29 


2122 


53. 


,13990 


-27.80500 


0. 


.180 


31.07 


0.91 


30.18 


0.28 


28.90 


0. 


.17 


29.32 


0.24 


2155 


53. 


,17774 


-27.80466 


0. 


.225 






29.14 


0.15 


27.49 


0. 


.06 


27.53 


0.07 


2276 


53. 


,16717 


-27.80429 


0. 


.126 






30.39 


0.34 


28.36 


0. 


.09 


28.21 


0.09 


2285 


53. 


,16834 


-27.80413 


0. 


,105 






29.54 


0.15 


27.71 


0. 


,04 


27.72 


0.05 


2303 


53. 


,18131 


-27.80394 


0. 


,159 






28.91 


0.11 


27.41 


0. 


.04 


27.44 


0.06 


2393 


53. 


,13875 


-27.80375 


0. 


,128 






34.19 


8.45 


29.12 


0. 


.13 


29.02 


0.14 
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Table 8 
First 36 I775-DROPOUT sources 



STScI ID RA (J2000) Dec (J2000) r50(") B cr^ V cr^ i cr^ z cr^ 



229 


53. 


17161 


-27, 


82079 


Q 


.160 






31 83 


]^ 


.47 


29 


04 





.13 


26, 


94 





.04 


263 


53 


17585 


-27 


81997 


Q 


239 


32 05 


4 24 








29 


47 





30 


27 


51 





12 


327 


53. 


16271 


-27. 


.81897 


0. 


,120 






32.82 


2 


.63 


29 


.54 


0. 


.14 


28, 


.14 





09 


344 


53. 


16691 


-27. 


.81875 


0. 


.066 


31.99 


1.34 


31.24 





.46 


29. 


90 


0. 


.15 


28. 


.13 


0, 


.07 


534 


53. 


16040 


-27. 


.81606 


Q 


230 


31.99 


3.35 


33.36 


g 


.11 


29 


.76 





.33 


27 


.30 





09 


591 


53. 


15532 


-27. 


.81516 


Q 


.150 


33.20 


7.34 


31.18 





.79 


32 


.08 


2 


.03 


27 


.31 





.07 


777 


53 


17048 


-27 


81243 


Q 


124 


30 84 


65 








29 


14 







27 


75 





08 


1148 


53. 


15203 


-27. 


80899 


0. 


.156 






31.42 





.85 


29 


32 


0. 


.14 


27, 


.89 


0, 


.10 


1149 


53. 


17985 


-27, 


80821 


Q 


.112 






32 14 




.28 


29 


91 





.19 


28, 


38 







1283 


53. 


14169 


-27 




Q 


,173 












29 


29 





.16 


27, 


85 





.11 


1309 


53. 


14225 


-27 


80677 


Q 


.098 


31.65 


1.15 








30 


.45 





30 


28, 


43 





10 


1312 


53 


13888 


-27 


80671 


Q 


108 


33 67 


8 61 








30 


53 





37 


28 


31 





14 


1526 


53. 


16672 


-27, 


.80416 


0. 


.157 






29.63 





.33 


26, 


.68 


0. 


02 


25, 


.11 


0, 


01 


1690 


53. 


16574 


-27, 


.80338 


0. 


.128 


32.39 


2.81 


31.13 


0, 


.61 


29. 


.24 


0. 


.12 


27. 


,84 


0, 


.08 


1790 


53. 


13043 


-27, 


.80231 


0. 


.167 


30.83 


0.91 


31.81 


1, 


.53 


28. 


,90 


0. 


.12 


27. 


,26 


0, 


.07 


1793 


53. 


14623 


-27, 


.80253 


0. 


,112 












30, 


.02 


0, 


.18 


28, 


.60 


0. 


,12 


1813 


53. 


,17748 


-27, 


.80245 


0. 


,126 












29, 


.56 


0, 


.17 


27, 


.62 


0. 


,08 


1852 


53. 


,14075 


-27, 


.80211 


0. 


.098 






33.24 


3 


.95 


29, 


.20 


0. 


.11 


27, 


.31 


0, 


.04 


1874 


53. 


,14617 


-27, 


.80190 


0. 


,137 






30.90 





.53 


30 


.01 


0. 


,26 


27, 


.91 


0, 


.10 


2007 


53. 


,17832 


-27, 


.80090 


0. 


.106 












31, 


.94 


1, 


.27 


28, 


.02 


0. 


.08 


2279 


53. 


14394 


-27, 


.79889 


0. 


.118 






31.84 


1 


.34 


28, 


.64 


0. 


.08 


27, 


.00 


0, 


.04 


2330 


53. 


,13587 


-27, 


.79833 


0. 


.145 


31.51 


1.54 


31.2 





.81 


28. 


.13 


0. 


.05 


26. 


.72 


0, 


.03 


2386 


53. 


14294 


-27, 


.79821 


0. 


.166 






30.95 





.59 


29, 


.16 


0. 


.13 


27, 


.71 


0, 


.08 


2408 


53. 


,14287 


-27, 


.79800 


0. 


.223 












29, 


.05 


0, 


.16 


27, 


.24 


0. 


.08 


2413 


53. 


,14285 


-27, 


.79785 


0. 


,192 


32.14 


3.01 


32.78 


3 


.81 


28. 


.63 


0. 


,09 


27. 


,28 


0, 


.06 


2490 


53. 


15950 


-27, 


.79758 


0. 


.127 






32.28 


1 


.92 


29, 


.07 


0. 


.11 


27, 


.46 


0, 


.06 


2521 


53. 


15261 


-27, 


.79724 


0. 


.314 












29, 


.38 


0, 


.27 


26, 


.99 


0. 


.08 


2760 


53. 


13928 


-27, 


.79580 


0. 


.112 


31.25 


0.97 


31.33 


0, 


.74 


29. 


.48 


0. 


.15 


27. 


,38 


0, 


.05 


2789 


53. 


17322 


-27, 


.79562 


0. 


,162 












29, 


.52 


0, 


.20 


27, 


.01 


0. 


.05 


2966 


53. 


,12271 


-27, 


.79459 


0. 


.200 


30.65 


0.79 


30.5 





.47 


29. 


.13 


0. 


.15 


27. 


,22 


0, 


.08 


3003 


53. 


14607 


-27, 


.79449 


0. 


,120 












31, 


.43 


0, 


.85 


27, 


.81 


0. 


.08 


3527 


53. 


12559 


-27, 


.79124 


0. 


.236 






31.05 





.95 


29, 


.13 


0. 


.18 


27, 


.42 


0, 


.09 


4266 


53. 


16193 


-27, 


.78699 


0. 


.076 


31.09 


0.61 


31.80 


0, 


.82 


30. 


.09 


0. 


.18 


28. 


,50 


0, 


.11 


4321 


53. 


18625 


-27, 


.78655 


0. 


,185 


29.90 


0.38 


30.68 


0, 


.56 


29. 


.63 


0. 


.24 


27. 


,43 


0, 


.09 


4353 


53. 


16161 


-27, 


.78633 


0. 


,186 












29, 


.66 


0, 


.21 


27, 


.86 


0. 


.09 


4715 


53. 


,13482 


-27, 


.78412 


0. 


.180 


31.66 


1.85 








29, 


.84 


0. 


.27 


27, 


.73 


0, 


.11 



Table 9 

Hubble dropout selection volumes 



B435-dropouts Veoe-dropouts irTs-dropouts 



Spectrum 


Zmin 




V(Mpc^) 


Zmin 




V(Mpc^) 


Zmin 


ZHmax 




ZGmax 


Vg(Mpc3) 


10 Myr 


3.53 


4.80 


3981 


4.67 


5.71 


2990 


5.71 


7.29 


3924 


6.62 


2335 


100 Myr 


3.48 


4.74 


4046 


4.62 


5.68 


3051 


5.68 


7.24 


3897 


6.56 


2279 


1000 Myr 


3.49 


4.69 


3866 


4.60 


5.65 


3037 


5.65 


7.21 


3893 


6.52 


2253 


Flat Step 


3.66 


4.75 


3466 


4.72 


5.70 


2797 


5.70 


7.19 


3734 


6.47 


2007 


Red Step 


3.72 


4.59 


2787 


4.69 


5.64 


2717 


5.64 


7.14 


3785 


6.42 


2056 




V is 


the co-moving volume (Mpc*^ arcmin ^) 


within the redshift range 


: V = 


arcmin^ r^max 







Vg and Vh are GOODS and HUDF, respectively. ZHmax ZGmax illustrate — — 20.8 at the zgso limits. 
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Table 10 

Hubble dropout selection volumes: Eb-v = 0.15 



B435-dropouts Veoe-dropouts irrs-dropouts 



Spectrum 




Zmax 


V(Mpc^) 




Zmax 


V(Mpc^) 




ZHmax 


V^f(Mpc3) 


ZGmax 


Vg(Mpc3) 


10 Myr 


3.35 


4.46 


3614 


4.53 


5.58 


3046 


5.58 


7.17 


4031 


6.58 


2456 


100 Myr 


3.30 


4.40 


3622 


4.50 


5.55 


3061 


5.55 


7.13 


4011 


6.55 


2404 


1000 Myr 


3.31 


4.33 


3385 


4.48 


5.52 


3031 


5.52 


7.10 


4022 


6.43 


2398 


Flat Step 


3.46 


4.36 


2956 


4.56 


5.55 


2869 


5.55 


7.08 


3870 


6.37 


2161 


Red Step 


3.51 


4.05 


1804 


4.55 


5.48 


2719 


5.48 


7.03 


3972 


6.33 


2258 




V is 


the co-moving volume (Mpc*^ arcmin ^) 


within the redshift range 


: V = 


arcmin^ r^^macc 







Vg and Vh are GOODS and HUDF, respectively. ZHmax & zcmax illustrate = —20.8 at the zgso limits. 



Table 11 

Surface densities of B435-DROPOUTS vs. magnitude 



Veoe range 


S(HUDF) 


S(GOODS-N) 
mag"""^ arcmin"^ 


S(GOODS-S) 


24-25 




0.069d=0.024 


0.094±0.029 


25 26 




0.808±0.155 


1.021±0.191 


26 26.5 




2.627±0.482 


3.031±0.551 


26.5 26.75 




3.966±0.745 


4.585±0.85 


26.75 27 




5.633±1.029 


5.286±0.97 


27 27.5 




6.479±1.138 


6.188±1.088 


25 26 


0.727±0.285 






26 26.5 


4. ±1.091 






26.5 27 


4.727±1.227 






27 27.5 


6.545±1.558 






27.5-28 


12.364±2.582 






28-28.5 


19.273±3.773 






28.5 29 


19.818±3.867 






29-29.3 


25.758±5.194 






29.3-29.624 


12.065±2.755 







Total surface brightness derived from HUDF only: 
23.7 mag arcmin"^, at Arest ~ 1400 A. 



Table 12 

Surface densities of V606-dROpouts vs. magnitude 



i775 range 


S(HUDF) 


S(GOODS-N) 
mag~^ arcmin"^ 


S(GOODS-S) 


24-25 




0.025±0.013 


0.025±0.013 


25 26 




0.295±0.066 


0.319±0.07 


26-26.5 




1.091±0.219 


1.006±0.204 


26.5-26.75 




1.731±0.361 


1.374±0.297 


26.75-27 




2.083±0.421 


1.202±0.267 


27-27.5 




2.371±0.438 


1.472±0.284 


25-26 


0.455±0.217 






26-27 


0.727±0.285 






27-27.5 


4. ±1.091 






27.5-28 


4.182±1.125 






28-28.5 


9.091±2.01 






28.5-29 


10.182±2.202 






29-29.25 


9.818±2.521 






29.25 29.524 


4.313±1.403 







Total surface brightness derived from HUDF only: 
24.8 mag arcmin"^, at Arest ~ 1400 A. 
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Table 13 

Surface densities of 1775-DROPOUTS vs. magnitude 



Z850 range 


S(HUDF) 


S(GOODS-N) 
mag""'^ arcmin"^ 


S(GOODS-S) 


25 26 




0.1±0.03 


0.081±0.026 


26 26.5 




0.339±0.087 


0.401±0.098 


26.5 27 




0.84±0.176 


0.915±0.189 


25 26.7 


0.214±0.113 






26.7 27.7 


2.d=0.545 






27.7-28.746 


2.434±0.619 







Total surface brightness derived from HUDF only: 
25.6 mag arcmin"^, at Arest ~ 1400 A. 



Table 14 
Schechter function fits 



Sample 


m* 


Simple Schechter 
S,(min-2) 




, (f)^ from 
M* (equiv) c 


(equiv) 


M* 


(/)*(Mpc~^) 




= -20.3 

(/>*(Mpc-3) x' 


B435-dropouts 


26.0 


6.6 


0.5 


-20.5 


0.0016 


-20.3 


0.0020 


0.4 


0.0020 0.4 


Veoe-dropouts 


26.1 


2.4 


1.2 


-20.6 


0.0008 


-20.5 


0.0010 


1.1 


0.0012 1.3 


i775-dropouts 


26.3 


1.5 


0.8 


-21.3 


0.0004 


-20.5 


0.0008 


0.9 


0.0010 1.1 



Columns 5 & 6 translate the simple Schechter function fits from columns 2 and 3 into equivalent and 0* assuming 
spectra from 100 My r CSF galaxy as in other cases. All contours are highly elliptical in the m*, Yj^ or M^,cf)^ planes, 
indicating strong degeneracy in the fitting parameters. All fits assume a = —1.6. The uncertainties in m^, are typically: 



cr(B435-drop): A(m(M), - ln[S,(0,)] ^ 0.6, A(m(M), + ln[S, (0, )] ^ 0.1; 
cr(V606-drop): A(m(M), - Inp, ((/>/)] ^ 1.4, A(m(M), + ln[S, ((/>/)] ^ 0.3; 
cr(i775-drop): A(m(M)* - ln[Y ^ {(f) ^)] ^ 1.5, A(m(M)\, + Inp* (0*)] ^ 0.5 



Table 15 

Combined luminosity functions G- and B435-DROPOUTS 



V 


mag 




log((/)) (mag' 


Mpc-^) 










SAGDP99 


HUDF 


GOODS-N 


GOODS-S 


-22.95 


to 


-22.45 


-5.48±0.18 








-22.45 


to 


-21.95 


-4.93±0.15 








-21.95 


to 


-21.45 


-4.24d=0.08 








-21.45 


to 


-20.95 


-3.72±0.05 








-22.48 


to 


-21.48 






-4.72±0.13 


-4.58±0.12 


-21.48 


to 


-20.48 




-3.71±0.14 


-3.63±0.08 


-3.52±0.07 


-20.48 


to 


-19.98 




-2.96±0.10 


-3.10±0.07 


-3.03±0.07 


-19.98 


to 


-19.73 






-2.92±0.07 


-2.85±0.07 


-19.73 


to 


-19.48 






-2.75±0.07 


-2.77±0.07 


-19.98 


to 


-19.48 




-2.88±0.10 






-19.48 


to 


-18.98 




-2.73d=0.09 


-2.60±0.07 


-2.62±0.07 


-18.98 


to 


-18.48 




-2.44±0.08 






-18.48 


to 


-17.98 




-2.23±0.08 






-17.98 


to 


-17.48 




-2.23±0.08 






-17.48 


to 


-17.18 




-2.10±0.08 






-17.18 


to 


-16.86 




-2.32±0.09 







Sample redshifts: 4.13 (3.75-4.5) Caltech, 4.11 (3.49-4.69) Hubble 
Sample volumes (Mpc^ arcmin"^): 1999 (Caltech), 4045 (Hubble) 
100 Myr Bruzual model distance moduU: 45.95 Caltech, 46.48 Hubble 
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Table 16 
Luminosity Density Comparison 









5](mag 


arcmin ^) 


log <l>(erg s ^ 


Hz-^ Mpc-^) 


^fair^t /^bright 


Sample 


z 




GOODS m < 


HUDF m > 


GOODS m < 


HUDF m > 




B435- dropout 


4.11 


25.31 


26.53±0.12 


24.49±0.05 


25.18±0.05 


26. ±0.02 


6.6±13% 


Veoe-dropout 


5.04 


25.66 


27.25±0.15 


25.57±0.07 


25.08±0.06 


25.75±0.03 


4.7±16% 


i775-dropout 


6.46 


26.06 


28.15±0.18 


25.88±0.14 


24.85±0.07 


25.76±0.06 


8.1±22% 



Faint/bright divide at M = —20.8. Faint e nd limit at M = —18.8. Mean redshifts, volumes from the 100 Myr model in Table 9. 
Listed uncertainties are oc 1/ Vnobj — 1 and do not include cosmic error or model-dependent variations. 
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